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A Mathematical Theory of Communication 

By C. E. SHANNON 
Introduction 

THE recent development of various methods of modulation such as PCM and PPM which exchange 
bandwidth for signal-to-noise ratio has intensified the interest in a general theory of communication A 
basis for such a theory is contained in the important papers of Nyquist 1 and Hartley- on this subject. In the 
present paper we will extend the theory to include a number of new factors, in particular the effect of no.se 
in the channel, and the savings possible due to the statistical structure of the original message and due to the 
nature of the final destination of the information. .... 

The fundamental problem of communication is that of reproducing at one point either exactly or ap- 
proximately a message selected at another point. Frequently the messages have meaning; that is they refer 
to or are correlated according to some system with certain physical or conceptual entities. These semantic 
aspects of communication are irrelevant to the engineering problem. The significant aspect is that the actual 
message is one selected from a set of possible messages. The system must be designed to operate for each 
possible selection, not just the one which will actually be chosen since this is unknown at the time of design. 

If the number of messages in the set is finite then this number or any monotonic function of this number 
can be regarded as a measure of the information produced when one message is chosen from the set, all 
choices being equally likely. As was pointed out by Hartley the most natural choice is the logarithmic 
function Although this definition must be generalized considerably when we consider the influence of the 
statistics of the message and when we have a continuous range of messages, we will in all cases use an 
essentially logarithmic measure. 

The logarithmic measure is more convenient for various reasons: 

1 It is practically more useful. Parameters of engineering importance such as time, bandwidth, number 
of relays etc., tend to vary linearly with the logarithm of the number of possibilities. For example, 
adding one relay to a group doubles the number of possible states of the relays. It adds 1 to the base 2 
logarithm of this number. Doubling the time roughly squares the number of possible messages, or 
doubles the logarithm, etc. 

2 It is nearer to our intuitive feeling as to the proper measure. This is closely related to (1) since we in- 
tuitively measures entities by linear comparison with common standards. One feels, for example, that 
two punched cards should have twice the capacity of one for information storage, and two identical 
channels twice the capacity of one for transmitting information. 

3. It is mathematically more suitable. Many of the limiting operations are simple in terms of the loga- 
rithm but would require clumsy restatement in terms of the number of possibilities. 

The choice of a logarithmic base corresponds to the choice of a unit for measuring information. If the 
base 2 is used the resulting units may be called binary digits, or more briefly bits, a word suggested by 
J W Tukey. A device with two stable positions, such as a relay or a flip-flop circuit, can store onebit of 
information. N such devices can store N bits, since the total number of possible states is 2 and log 2 2 - N. 
If the base 10 is used the units may be called decimal digits. Since 

log 2 M = log 10 M/log 10 2 
= 3.321og 10 M, 

'Nyquist, H., "Certain Factors Affecting Telegraph Speed" Bell System Technical Journal, April 1924, p. 324; "Certain Topics in 
Telegraph Ti^missionTheory;M./.£.£. Trow., v. 47, April 1928, p. 617. 

'Hartley, R. V. L., 'Transmission of Information," Bell System Technical Journal, July 1928, p. 533. 
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Fig. 1 -Schematic diagram of a general communication system. 

consists of essentially five parts: 

1. An,/^^ 

receiving terminal. The message may be o « ype, (a) A s^u a rf 

of teletype system; (b) A r^SS^^-^^^^ 1 ^^ 
time and other variables as in black and white television ( on a 

function /(*,y,0 of two space coordinates and time the ligh mtens.ry atpo^ , W fa ^ 
pickup tube plate; (d) Two or more functions of T?"^'*^]^ individual channels in 
"dimensional" sound transmission or if the »^ B ^ to tt ^^ lW co 0 «*rf^ 
nmlt.plex;^^ 

functions /(x.y.r). *(*.*') d f fi T denn ed in the region - similarly, several 

of these three functions as components of a vector , t ^™~ 1 * f a number of functions 

black and white television sources would I«f«JJ^^ Virion with an associated 
of three variables; (f) Various combinations also occur, for example 

audio channel. 

2. A — which operates on the me*.* SS—^ 
mission over the channel. In telephony h « opmdort cowtt me «y 8 ^ 
into a proportional electrical cunent. In <**^J™ZZ?* me message. In a multiplex 
. sequence of dots, dashes and spaces on the channel ^y™"*n S sed, and encodedi 

PCM system me different speech f"nc« '^^^L^m and frequency 

the message from the signal. 
5 The destination is the person (or thing) for whom the message is intended. 
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physical counterparts. We may roughly classify communication systems mto three main categon^ d.^ 
continuous and mixed. By a discrete system we will mean one in which both the message and the ^gnal 
are a sequence of discrete symbols. A typical case is telegraphy where the message ,s a sequence o letters 
and the signal a sequence of dots, dashes and spaces. A continuous system ,s one in which the message and 
signal are both treated as continuous functions, e.g., radio or television. A m.xed system .s one in winch 
both discrete and continuous variables appear, e.g., PCM transmission of speech. _ 

We first consider the discrete case. This case has applications not only in communicanon theory but 
also in the theory of computing machines, the design of telephone exchanges and other fields. In addition 
the discrete case forms a foundation for the continuous and mixed cases which will be treated m the second 
half of the paper. 

PART I: DISCRETE NOISELESS SYSTEMS 
1. The Discrete Noiseless Channel 

Teletype and telegraphy are two simple examples of a discrete channel for transmitting information. Gen- 
m lya d screte channel will mean a system whereby a sequence of choices from a finite set of elementary 
symbols S, , . . . , S n can be transmitted from one point to another. Each of the symbols S,- .s assume^ tohave 
a certain duration in time f, seconds (not necessarily the same for different S h for example the dots and 
dashe in telegraphy). It is not required that all possible sequences of the ,* be capable of «,s sion on 
the system; certain sequences only may be allowed. These will be possible signals for he channeL Thus 
in telegraphy suppose the symbols are: (1) A dot, consisting of line closure for a umt of tune and then line 
open for aVnit oftime; (2) A dash, consisting of three time units of closure and one unit open; (3) A letter 
space consisting of, say, three units of line open; (4) A word space of six units of line open. We 
I restriction on allowable sequences that no spaces follow each other (for if two letter spaces are adjacent 
it is identical with a word space). The question we now consider is how one can measure the capacty of 

such a channel to transmit information. _ f , K „ v> c„mhnk 

In the teletype case where all symbols are of the same duration, and any sequence of the 32 symbols 
is allowed the answer is easy. Each symbol represents five bits of information. If the system transmits n 
symbols per second it is natural to say that the channel has a capacity of 5n bits per second. This does not 
mean that the teletype channel will always be transmitting information at this rate - this is the maximum 
possible rate and whether or not the actual rate reaches this maximum depends on the source of mformation 
which feeds the channel, as will appear later. 

In the more general case with different lengths of symbols and constraints on the allowed sequences, we 

make the following definition: 

Definition: The capacity C of a discrete channel is given by 

where N(T) is the number of allowed signals of duration T. >u ot tht > t imit 

It is easily seen that in the teletype case this reduces to the previous result. It can be shown that the hm 
in question will exist as a finite number in most cases of interest. Suppose all sequences of the symbo 
S h ...,S„ are allowed and these symbols have durations t u ...,t n . What 1S the channel capacity? If N(t) 
represents the number of sequences of duration t we have 

N(t) =N{t-ti)+N{t-t 2 ) + --+ N(t - 1„). 

The total number is equal to the sum of the numbers of sequences ending in S, S 2 , ... , J .and i ftese > are 
N(t - r, ) Nit - ti) N(t - 1„), respectively. According to a well-known result in finite differences, N{t) 
is then asymptotic for large t to where Xo is the largest real solution of the characteristic equation: 

x -r 1+x -«2 + ... + x-'" = l 
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and therefore C = logX 0 . 

type and find C from the characteristic equation. In the telegraphy 

w(f) «^-2)+Nfr-4)+N(»-5)+N(.-7) + N(.-D+N(#-lfl) 

as we see by counting sequences of symbols according ; to .the f^^J^Z^ 
HenceCis-lbg^wherewistheposmverootof^M +V +/* + ^ 

C = 0.539. . . .„ , _ n a ii OW ed sequences is the following: We 

A very general type of restriction which may b ^^^^^^^^Su^ 
imagineanumberofpossiblestatesa,,^ ^^.^T^Zom of these has been transmitted the 
can be transmitted (different the particular symbol transmitted. The 

state changes to a new state *P?^*^^^^ depending on whether or not a space was 
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Fig 2-Graphical representation of the constraints on telegraph symbols, 
in accordance with the following result: . 



where Sn = 1 if i = 7 and is zero otherwise. 

For example, in the telegraph case (Fig. 2) the determinant is: 

I _i (vv- 2 +w-*) [ o. 

| (VV -3 + W-6) (W- 2 + W" 4 -l)| 

On expansion this leads to the equation given above for this case. 

2. THE DISCRETE SOURCE OF INFORMATION 



Wn no„ consider one inforrnanon »^^^™"™„^, The main point a. ba" is the 



of proper encoding of the information. In telegraphy, for example, the messages to be transmitted consist of 
sequences of letters. These sequences, however, are not completely random. In general, they form sentences 
and have the statistical structure of , say, English. The letter E occurs more frequently than Q, the sequence 
TH more frequently than XP, etc. The existence of this structure allows one to make a saving in time (or 
channel capacity) by properly encoding the message sequences into signal sequences. This is already done 
to a limited extent in telegraphy by using the shortest channel symbol, a dot, for the most common English 
letter E; while the infrequent letters, Q, X, Z are represented by longer sequences of dots and dashes. This 
idea is carried still further in certain commercial codes where common words and phrases are represented 
by four- or five-letter code groups with a considerable saving in average time. The standardized greeting 
and anniversary telegrams now in use extend this to the point of encoding a sentence or two into a relatively 
short sequence of numbers. 

We can think of a discrete source as generating the message, symbol by symbol. It will choose succes- 
sive symbols according to certain probabilities depending, in general, on preceding choices as well as the 
particular symbols in question. A physical system, or a mathematical model of a system which produces 
such a sequence of symbols governed by a set of probabilities, is known as a stochastic process. 3 We may 
consider a discrete source, therefore, to be represented by a stochastic process. Conversely, any stochastic 
process which produces a discrete sequence of symbols chosen from a finite set may be considered a discrete 
source. This will include such cases as: 

1 . Natural written languages such as English, German, Chinese. 

2. Continuous information sources that have been rendered discrete by some quantizing process. For 
example, the quantized speech from a PCM transmitter, or a quantized television signal. 

3. Mathematical cases where we merely define abstractly a stochastic process which generates a se- 
quence of symbols. The following are examples of this last type of source. 

(A) Suppose we have five letters A, B, C, D, E which are chosen each with probability .2, successive 
choices being independent. This would lead to a sequence of which the following is a typical 
example. 

BDCBCE C C CADCBDDA A E C E E A 
ABBDAEECACEEB AEEC BCEAD. 

This was constructed with the use of a table of random numbers. 4 

(B) Using the same five letters let the probabilities be .4, .1, .2, .2, .1, respectively, with successive 
choices independent. A typical message from this source is then: 
AAACDCBDCEAADADACEDA 
EADCABEDADDCECAAAAAD. 

(C) A more complicated structure is obtained if successive symbols are not chosen independently 
but their probabilities depend on preceding letters. In the simplest case of this type a choice 
depends only on the preceding letter and not on ones before that. The statistical structure can 
then be described by a set of transition probabilities pi{j), the probability that letter i is followed 
by letter The indices i and j range over all the possible symbols. A second equivalent way of 
specifying the structure is to give the "digram" probabilities p(i, i.e., the relative frequency of 
the digram i j. The letter frequencies (the probability of letter Q, the transition probabilities 

3 See, for example, S. Oiandrasekhar, "Stochastic Problems in Physics and Astronomy," Reviews of Modern Physics, v. 15, No. 1, 
January 1943, p. 1. 

4 Kendall and Smith, Tables of Random Sampling Numbers, Cambridge, 1939. 
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Pi (j) and the digram probabilities p(«, j) are related by the following formulas: 

p(0 = 2pM = 2>O'.0 = 5>0")pj(0 

7 i ' 

p(«,;) = p(')piU) 
Sp,O) = Xp(0 = Sp(^) = 1 - 

As a specific example suppose there are three letters A, B, C with the probability tables: 



PiU) 


ABC 


p(i) p(U) 


ABC 


A 
i B 
C 


0|! A 
I. i 0 B 

1 3 TO C 


A A 

1 * c 


0 £ 13 
A * 0 

* -4. i 

27 135 135 



A typical message from this source is the following: ab hhraCACAB 

ABBABABABABABABBBABBBBBABABABABABBBACACAB 

B ABBBB ABB AB ACBBB AB A. . 
The next increase in complexity would involve trigram frequencies but no more. The choiceof 

I ^ would depend on'the preceding two letters but not on 

set of trigram frequencies p(i,j,k) or equivalently a set of ^^J^TJ^^. 
be required. Continuing in this way one obtains « sa ^. f ? m ( ^^^^ 
cesses In the general /.-gram case a set of n-gram probabilities p{iun,-^) or of transition 

probabilities p nh ,„_, (<„) is required to specify the statistical structure. 

(D) Stochastic processes can also be definec .which produce a - — of „e of 
"words." Suppose there are five letters A, r>, U u, n ami 
associated probabilities: 

10 A 16 BEBE .11 CABED .04 DEB 

04ADEB .04 BED .05 CEED .15 DEED 
05ADEE .02 BEED .08 DAB .01 EAB 
.01 BADD .05 CA .04 DAD .05 EE 
Suppose successive "words" are chosen independently and are separated by a space. A typical 

SSS A mBE DEED DEB ADEE ADEE EE DEB BEBE BEBE BEBE ADEE BED DEED 
SS S dEED DEED BEBE CABED BEBE BED DAB DEED ADEB^ 
If all the words are of finite length this process is equivalent to one 
the description may be simpler in terms of the word structure and probability. We may 
generalize here and introduce transition probabilities between words, etc. 

each letter having the same probability that it has m the w wim 
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u to P n etc In the second-order approximation, digram structure is introduced. After a 

letters. 

3 The Series of approximations to English 

"alphabet," the 26 letters and a space. 

1 Zero-order approximation (symbols independent and equiprobable). 

XFOML RXKHRJFFJUJ ZLPWCFWKCYJ FFJEYVKCQSGHYD QPAAMKBZAACBZL- 
HJQD. 

2 First-order approximation (symbols independent but with frequencies of English text). 

OCRO HLI RGWR NMIELWIS EU LL NBNESEBYA TH EEI ALHENHTTPA OOBTTVA 
NAHBRL. 

3 Second-order approximation (digram structure as in English). 

ON IE ANTSOUTINYS ARE T INCTORE ST BE S DEAMY ACHIN D ILONASIVE TU- 
CoS^T TEASONARE FUSO TIZIN ANDY TOBE SEACE CTISBE. 

4 Third-order approximation (trigram structure as in English). 

IN NO 1ST LAT WHEY CRATICT FROURE BIRS GROCID PONDENOME OF DEMONS- 
TURES OF THE REPTAGIN IS REGOACTIONA OF CRE. 

appropriate frequencies. 

THE LINE MESSAGE HAD BE THESE. 
6. Second-order word approximation. The word transition probabilities are correct but no further struc- 
ture is included. 

THE^TME ^WHOEVER TOLD THE PROBLEM FOR AN UNEXPECTED. 
TheresemblancetoordinaryEnglishtextincreases^ 
these samples have reasonably good structure out to about m *e 3J£™ nces> but fou , 
construction. Thus in (3) the statistical process msur -"^^T^ sequences of four or more 
letter sequences from the sample can usually be fitted ^J^^^^ particular sequence 
words can easily be placed ,in sentences ^^^ a e ^SZ a JL.^ V ^^ 
of ten words "attack on an Enghsh wnter that the xhmcter of thu mi ™ discrete source . 



To construct (3) for example, one opens a book at random and selects a letter at random on the page. This 
letter is recorded. The book is then opened to another page and one reads until this letter is encountered. 
The succeeding letter is then recorded. Turning to another page this second letter is searched for and the 
succeeding letter recorded, etc. A similar process was used for (4), (5) and (6). It would be interesting if 
further approximations could be constructed, but the labor involved becomes enormous at the next stage. 

4. Graphical Representation of a Markoff Process 

Stochastic processes of the type described above are known mathematically as discrete Markoff processes 
and have been extensively studied in the literature. 6 The general case can be described as follows: There 
exist a finite number of possible "states" of a system; S U S 2 , - - • ,5„. In addition there is a set of transition 
probabilities; pifj) the probability that if the system is in state S; it will next go to state Sj. To make this 
Markoff process into an information source we need only assume that a letter is produced for each transition 
from one state to another. The states will correspond to the "residue of influence" from preceding letters. 
The situation can be represented graphically as shown in Figs. 3, 4 and 5. The "states" are the junction 




d .2 

Fig. 3 — A graph corresponding to the source in example B. 

points in the graph and the probabilities and letters produced for a transition are given beside the correspond- 
ing line. Figure 3 is for the example B in Section 2, while Fig. 4 corresponds to the example C. In Fig. 3 




Fig. 4 — A graph corresponding to the source in example C. 

there is only one state since successive letters are independent. In Fig. 4 there are as many states as letters. 
If a digram example were constructed there would be at most n 2 states corresponding to the possible pairs 
of letters preceding the one being chosen. Figure 5 is a graph for the case of word structure in example D. 
Here S corresponds to the "space" symbol. 

5 . Ergodic and Mixed Sources 

As we have indicated above a discrete source for our purposes can be considered to be represented by a 
Markoff process. Among the possible discrete Markoff processes there is a group with special properties 
of significance in communication theory. This special class consists of the "ergodic" processes and we 
shall call the corresponding sources ergodic sources. Although a rigorous definition of an ergodic process is 
somewhat involved, the general idea is simple. In an ergodic process every sequence produced by the process 

6 For a detailed treatment see M. Frechet, Mfthode des Sanctions arbilraires. Thforie des tenements en chaine dans le cos d'un 
nombre fini d'itats possibles. Paris, Gauthier-Villars, 1938. 



^JU »>. Roughly ft. ergodic property means s«H is relaed to *« smicu-e 

impossible to go from junctions in part B to junctions in part A. 

2 . A closed series of lines in the graph with all arrows — — EES 

be called a "circuit." The "length" of a circuit .s the '^^^^^^^ 
is a circuit of length 5. The second property required ,s that the greatest common a 
of all circuits in the graph be one. 




Fig. 5^A graph corresponding to the source in example D. 

Ifthenrstconditionissati^butthesecondo^ 
to d > 1, the sequences have a certain type of penod. The .™ s ^ , .„ ^ e is 

classes which are sta^^^ 

called letter 1) By a shift of from 0 up to d J W «*« ^ a b c Letter a is 

a Thus a typical sequence is 

This type of situation is not of much importance for our work. . achofwhichsatisfies 

If te first condition is violated the graph may be '^^^^^ wi ^^^ 
thenrstcondi.cn. We^u^ 

'THese are restatements in terms of the graph of conditions given in Frechet 
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where p\ is the probability of the component source L,\ 

Physically the situation represented is this: There are several different sources L\ 9 L 2 , Li, . . - which are 
each of homogeneous statistical structure (i.e., they are ergodic). We do not know a priori which is to be 
used, but once the sequence starts in a given pure component L,-, it continues indefinitely according to the 
statistical structure of that component. 

As an example one may take two of the processes defined above and assume p\ = .2 and p 2 = .8. A 
sequence from the mixed source 

L=.2Li + .8L 2 

would be obtained by choosing first L\ or L 2 with probabilities .2 and .8 and after this choice generating a 
sequence from whichever was chosen. 

Except when the contrary is stated we shall assume a source to be ergodic. This assumption enables one 
to identify averages along a sequence with averages over the ensemble of possible sequences (the probability 
of a discrepancy being zero). For example the relative frequency of the letter A in a particular infinite 
sequence will be, with probability one, equal to its relative frequency in the ensemble of sequences. 

If Pi is the probability of state i and pi(j) the transition probability to state j, then for the process to be 
stationary it is clear that the Pi must satisfy equilibrium conditions: 

i 

In the ergodic case it can be shown that with any starting conditions the probabilities Pj{N) of being in state 
j after N symbols, approach the equilibrium values as N -> °°. 

6. Choice, Uncertainty and Entropy 

We have represented a discrete information source as a Markoff process. Can we define a quantity which 
will measure, in some sense, how much information is "produced" by such a process, or better, at what rate 
information is produced? 

Suppose we have a set of possible events whose probabilities of occurrence are pi,p 2 ,. • - ,Pn- These 
probabilities are known but that is all we know concerning which event will occur. Can we find a measure 
of how much "choice" is involved in the selection of the event or of how uncertain we are of the outcome? 

If there is such a measure, say H{p\ ,p 2 , . • ■ ,P«)> it is reasonable to require of it the following properties: 

1 . H should be continuous in the p,-. 

2. If all the p\ are equal, p, = then H should be a monotonic increasing function of n. With equally 
likely events there is more choice, or uncertainty, when there are more possible events. 

3. If a choice be broken down into two successive choices, the original H should be the weighted sum 
of the individual values of H. The meaning of this is illustrated in Fig. 6. At the left we have three 




Fig. 6— Decomposition of a choice from three possibilities. 

possibilities p\ = 3, pi = 5, P3 = On the right we first choose between two possibilities each with 
probability \, and if the second occurs make another choice with probabilities \ , 3. The final results 
have the same probabilities as before. We require, in this special case, that 

H(i.i.i) = »(i.i)+i»(§.i)- 

The coefficient i is because this second choice only occurs half the time. 
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In Appendix 2, the following result is established: 

Theorem 2: The only H satisfying the three above assumptions is of the form: 

n 

H = -K^ Pi logpi 
i=l 

where K is a positive constant. 

This theorem, and the assumptions required for its proof, are in no way necessary for the present theory. 
It is given chiefly to lend a certain plausibility to some of our later definitions. The real justification of these 
definitions, however, will reside in their implications. 

Quantities of the form H = -£ p, \og Pi (the constant K merely amounts to a choice of a unit of measure) 
play a central role in information theory as measures of information, choice and uncertainty. The form of H 
will be recognized as that of entropy as defined in certain formulations of statistical mechanics where Pi is 
the probability of a system being in cell i of its phase space. H is then, for example, the H in Boltzmann s 
famous H theorem. We shall call H = - 5>; logp, the entropy of the set of probabilities p, , . . . , p„. If * is a 
chance variable we will write H(x) for its entropy; thus x is not an argument of a function but a label for a 
number, to differentiate it from H(y) say, the entropy of the chance variable y. 

The entropy in the case of two possibilities with probabilities p and q = 1 - p, namely 



H = -{p\ogp + q\ogq) 



is plotted in Fig. 7 as a function of p. 




' 0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 

P 

Fig. 7— Entropy in the case of two possibilities with probabilities p and ( 1 - p). 



The quantity H has a number of interesting properties which further substantiate it as a reasonable 
measure of choice or information. 

1. H = 0 if and only if all the p, but one are zero, this one having the value unity. Thus only when we 
are certain of the outcome does //vanish. Otherwise Wis positive. 

2. For a given n, H is a maximum and equal to logn when all the Pi are equal (i.e., -). This is also 
intuitively the most uncertain situation. 

8 See, for example, R. C. Tolman, Principles of Statistical Mechanics, Oxford, Clarendon, 1938. 
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j° inteVentiS H( X ,y) = -%p(iJ)^p{iJ) 



while 



H(x) = -£pMlogI>M 

U J 

It is easily shown that H(x,y) <H{x) + H(y) 

■ a „„„,wrie o(i i) = p(i)p(j))- The uncertainty of a joint event is 
with equality only if the events are independent (i.e., p(i,J) PWPUV 

less *an or equal to tine sum o^^ increases //. Thus if p, < P2 and 

, h ,11 „ • > 0 then H increases (except in the special case where this transfer- 

M;) "X;P(U)' 

Wedennethe — ^af^)-*.^^^-' te -^ rf *^ 
according to the probability of getting that particular x. That is 

H x (y) = -J,P(U)^PiU)- 

Pi (j) we obtain 

Hx (y) = -£p(i,;)logp('j) + 
= H{x,y)-H{x) 

or H(x,y) = HW+H,(y). 

Th euncertainty(o^^^ 
known. 

6. From 3 and 5 we have ^ + ^ ^ ^ ^ = ^ + 

Hence HWZHtfy). 

Tneuncertaintyofyisnev^^ 
events, in which case it is not changed. 
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7. The Entropy of an Information Source 
Consider a discrete source of the finite state type considered above. For each possible state , there will be a 
set of probabilities *(/) of producing the various possible symbols j. Thus there is an en*opy H> for each 
state. The entropy of the source will be defined as the average of these H, weighted in accordance with the 
probability of occurrence of the states in question: 

H = 

i 

This is the entropy of the source per symbol of text. If the Markoff process is proceeding at a definite time 
rate there is also an entropy per second 

i 

where ft is the average frequency (occurrences per second) of state i. Clearly 

H' = mH 

where m is the average number of symbols produced per second. H or H' measures the amount of informa- 
tion generated by the source per symbol or per second. If the logarithmic base is 2, they w.ll represent bus 

^ IfTuclTivTsrmbolt are independent then H is simply - I W log W where Pi is the probability of sym- 
bol . Suppose in Lease we consider a long message of N symbols. It will contain with high ^probabdrty 
Tut recurrences of the first symbol, p 2 N occurrences of the second, etc. Hence the probability of tins 
particular message will be roughly ^ ^ 

P — Pi Pi m "Pn 



or 



logp = W^p,logpi 

i 

log p=-NH 
logl/£ 

H is thus approximately the logarithm of the reciprocal probability of a typical long sequence divided by the 
number of symbols in ihe sequence. The same result holds for any source. Stated more precisely we have 
(see Appendix 3): 

Theorem 3: Given any e > 0 and S > 0, we can find an JV 0 such that the sequences of any length N > N 0 
fall into two classes: 

1. A set whose total probability is less than e. 

2. The remainder, all of whose members have probabilities satisfying the inequality 



logp 1 



N 



-H 



<S. 



In other words we are almost certain to have l^J. very close to H when N is large. 

A closely related result deals with the number of sequences of various probabilities Consider again the 
sequences of length N and let them be arranged in order of decreasing probability. We define »(«)»>* 
dte number we must take from this set starting with the most probable one ,n order to accumulate a total 
probability q for those taken. 
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Theorem 4: \ogn{q) „ 

Lim — - — = H 



2_HA '- u .« t tj anH H - can be determined by limiting operations directly from 

* srss^~ - - - — — 

G„ = -i£p(ft)io8<>(ft) 

p(B h Sj)/p{Bi) be the conditional probability of Sj after Bi. Let 

F N = -^p(Bi,Sj)\ogp Bi (Sj) 

*re me aum * oue, „- 6,0* B, of* - 1 .»d - - *>■ «" * * ' —** 

decreasing function of N, 

F N =NG N -(N-l)G N -\, 
1 N 

and Lim//-»po Fn = H. *«« m vim»rirwK to H can be obtained 

Theseresults are derived in Appendix S. They show that a series of F „ is the 

by considering only the statistical structure of Ihe source of the type 

better approximation. In fact F N is ^.«JX^i^S^ than N ** '« * ** 

discussed above. If there are no statistical influences extending o angedby a knowledge of 

The nuio of .he entropy of a source to the max mum veto « couUto e ^ ^ 

symbols will becallcd its rctoivc crropy. This ,s <««• ""f^a^Ltac, ot ord i„a^ English, 
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this neighborhood. One is by calculation of the entropy of the approximations to English. A secon Imrfhod 
S o Ste I c^in fraction of the letters from a sample of English text and then let someone attempt to 
rln 1S1 If *ey can be restored when 50% are deleted the redundancy must be greater than 50%. A 
third method depends on certain known results in cryptography. . , 

Two extremes of redundancy in English prose are represented by Basic Enghsh and by James Joyce s 
The Basic English vocabulary is limited to 850 words and the redundancy ,s very 
gh 2 expansion** occurs when a passage is translated into Bas.c Enghs^ Joyce 

on the othe hand enlarges the vocabulary and is alleged to achieve a compression of semantic content. 

CmTCw !■ ^ ated to * e existenCe ° f CrOSSW ° rd PUZZleS " ^ "7 n « 
™^nv ^enc^ of etters ts a reasonable text in the language and any two-dimensional array of letters 

the redundancy is too high the language imposes too many constants for large 
forms a crossword P^™f ana , is shows that if we assume the constraints .mposed 

nrCag " . J?. , S£!Si3S« naLe, large crossword puzzles are just possible when 
^50%. If the redundancy is 33%, three-dimensional crossword puzzles should be possible, 

etc. 

8 REPRESENTATION OF THE ENCODING AND DECODING OPERATIONS 
>wmt ; a second function* these two quantities. Thus a transducer can be described by two functions. 

y n = f(Xn,OLn) 
a n +\ =g{Xn,Ot n ) 

where 

x n is the input symbol, 

a n is the state of the transducer when the input symbol is introduced, 

y n is the output symbol (or sequence of output symbols) produced when x n is introduced if the state is a n . 
If the output symbols of one transducer can be identified with the input symbols of a second they can be 

the second will be called its inverse. 

TWem 7- The output of a finite state transducer driven by a finite state statistic^ source ,s a .tote 
source, Sft entropy (per unit time) less than or equa, to mat of the mput. If the transducer 
is non-singular they are equal. 

Letarepresent the state of the source, which produces a sequence of symbols and let / 3bethes ^ 
the^rwhichproducesanitsoutpu^^ 

u *u «^A»m ct^tP Qn^ce" of Dairs fa BY Two points in the space (a\J\) and (a2,P2), are conneciea oy 
Ce if afcan " uce an * wE changes A toV and this line is given the probability of *at * m this 
a e The ine is labeled with the block of yj symbols produced by the transducer. The entropy o fe. outpu 
™ h ^calculated as the weighted sum over the states. If we sum first on 0 each resulting term is less than or 
^S2™«dJ^to a, hence the entropy is not increased. If the transducer is non-smgular 

L first and second transducers respectively, then H[ > H' 2 > H>, = H[ and therefore //, - W 2 . 
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*"*' „<•> = ?lw-* 



where <jj> is the duration of the S * symbol leading from state i to state / and the B, satisfy 



/ (5) 

e^X B i w y 



mized at the channel capacity. 

9 THE FUNDAMENTAL THEOREM FOR A NOISELESS CHANNEL 
We will now justify our interpretation of H as the rate of generating information by proving that H deter- 

e£- £ sym b olspersecondover the channel where e is arhi.anly small. It is not pos.hle to .ansm.t at 



H C 



an average rate greater than — . 

The convex pan of dye ton,, ,ha. £ canno. be exceeded, ra«y be proved by noting ft* .he endopy 

also this entropy cannot exceed the channel capacity, nence n _ 

" T£Z'£ of d,. to, „«. be jj- I in fwo J J£^^^££ 

one containing less than 2^ members ana me probability less than m- As /V increases 

small when T is large, if we choose 

r= (! + T 

thentherewillbeasu—^ 

N and T are sufficiently large (however small A) an es are represe nted by larger 
„ here „ is -I. The mean _ of emission ,„ sne^e symbofs per second wd, daen be ^..r than 

[ (1 -«i + 4p[u-<§ + AMM~' 
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AsNincreasesMand.approachzero^ 

and suppose their probabilities are 

Arrange the messages of length Nm order of de re^.ng prc ry ^ ^ ^ g> p ^ 

p, > P2 > P3 • • • > Pn ft = « IJ ^» P ^ ffor message * is obtained by expanding P s as a bmary 
We first encode into a binary system. The binary mm* > tor me g 
number. The expansion is carried out to m, places, where m s is the g 

log 2 j<m s <l + \og 2 -. 

^^^^^^^^^^^^ 
codes. From these inequalities we have 

1 ^ ^ 1 

2^<P*< 2mj -r 

SI, at least * large, and ** binary '*V>" s '™^Z?fZlTcote. If the ehannel sequences Me 
have ■ 1 v 



But, 



iS(i*>s»2>*<s£( ,+, *S> 



and therefore, 1 



As /V increases G„ approaches H the entropy of the source and H' approaches H. ^ ^ ^ ^ 
We see from this that the inefficiency in coding ™ on ^ and ^ entropy G * calculated for 



H HN 



This method of encoding is subs— y £ — . - ^ L^Tnlmo 
method is to arrange the messages of d ^ age f s P in & first group its first binary digit 

groups of as nearly equal probability as poss.We. * Qf nearly ^ probability and the 

will be 0, otherwise 1. The groups are singly J** „ continued unti , each subset contains 

oTyT^^ 

I the same thing as the arithmetic process described above. 

10 Discussion and Examples 

introduced so that the generator as seen from the toad has tn > Jn a staUstlcal 

»Tech»ica. Report No. 65, Ttie Research Laboratory of E.ecTonics, HIT, March .7, .949. 
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asthe sour ce ^^^^SS:^^ 
exact match is not in general poss.ble, we can W^JJ^wrf i coding system. This is of course 
rate of transmission to the capacity C may be called ^JJSto^Smimi possible entropy, 
equal totheratiooftheactualentropy o .thechann ^^^"an^itter and receiver. In the 
in general, ideal or nearly ideal encoding "J^J^^ ^ ^ is to allow reasonably good 
noiseless case which we have been constdenng * ma, ^nfuncuon ofth s y ^ ^ rf ^ 



fact j logp"' . c 



must be small for all but a small * *» is zer0 , and n0 channel is required. For 

If a source can produce only one particu ar me* a ge £"JJW duces a definite sequence 

example, a computing machine set up to ^^^^^ One could construct a 
with no chance element. No channel is required o — Q *^ ismay L impractical. In such a ca^e 
•econdm-cMne^ 

we can choose to ignore some or all of the » *^J£J£ J capab i e of sending any sequence of 

the digits of .r to be a random sequence in ^J^^^d know ledge of English in constructing 
digits In a similar way we may choose to use some > of o* ^ entropy subject t0 the 

a code, but not all of it. In such a case we confer ^^^^ ^ channel capacit y which 
statistical conditions we w.sh to retain. The enttopy 1 0 ™ d „ mat al , the digits are chosen 

is necessary and sufficient. In the , «^*^^^the statistical saving possible due to 

and is3e d =?^ 

cho S enfromamongA,B,C,D with probabilities 5 , 3 , s , 8 , successive sy 

We have ; 

H = -Qlog^ + ^ogJ + flog 8 ) 

= | bits per symbol. 

code (obtained by the method of the second proof of Theorem 9) . 

A 0 
B 10 

c no 

D HI 

The average number of binary digits used in encoding a sequence of N symbols will be 

W(Jxl + ix2+|x3) = JiV. 

bit per symbol. Since, on the average, we have , btn»y! is | og4 = 2, oocorring when n. B, C, 

b „t *- as. 1 ? nrsss. sr^TS. - b - ^ 10,0 

00 A' 

01 B' 

10 C 

11 D> 
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This double process then encodes the original message into the same symbols but with an average compres- 
sion ratio g. 

As a second example consider a source which produces a sequence of A's and ZTs with probability p for 
A and q for B. If p «: q we have 

H = -\ogp p {l-p) l - p 
= -p\ogp(\-pf-M p 

= plog-. 
P 

In such a case one can construct a fairly good coding of the message on a 0, 1 channel by sending a special 
sequence, say 0000, for the infrequent symbol A and then a sequence indicating the number of B's following 
it. This could be indicated by the binary representation with all numbers containing the special sequence 
deleted. All numbers up to 16 are represented as usual; 16 is represented by the next binary number after 16 
which does not contain four zeros, namely 17 = 10001, etc. 

It can be shown that as p -> 0 the coding approaches ideal provided the length of the special sequence is 
properly adjusted. 

PART II: THE DISCRETE CHANNEL WITH NOISE 

1 1 . Representation of a Noisy Discrete Channel 

We now consider the case where the signal is perturbed by noise during transmission or at one or the other 
of the terminals. This means that the received signal is not necessarily the same as that sent out by the 
transmitter. Two cases may be distinguished. If a particular transmitted signal always produces the same 
received signal, i.e., the received signal is a definite function of the transmitted signal, then the effect may be 
called distortion. If this function has an inverse — no two transmitted signals producing the same received 
signal — distortion may be corrected, at least in principle, by merely performing the inverse functional 
operation on the received signal. 

The case of interest here is that in which the signal does not always undergo the same change in trans- 
mission. In this case we may assume the received signal E to be a function of the transmitted signal S and a 
second variable, the noise N. 

E = f(S ) N) 

The noise is considered to be a chance variable just as the message was above. In general it may be repre- 
sented by a suitable stochastic process. The most general type of noisy discrete channel we shall consider 
is a generalization of the finite state noise-free channel described previously. We assume a finite number of 
states and a set of probabilities 

This is the probability, if the channel is in state a and symbol i is transmitted, that symbol j will be received 
and the channel left in state /3. Thus a and f3 range over the possible states, i over the possible transmitted 
signals and j over the possible received signals. In the case where successive symbols are independently per- 
turbed by the noise there is only one state, and the channel is described by the set of transition probabilities 
Pi(j), the probability of transmitted symbol i being received as 

If a noisy channel is fed by a source there are two statistical processes at work: the source and the noise. 
Thus there are a number of entropies that can be calculated. First there is the entropy H{x) of the source 
or of the input to the channel (these will be equal if the transmitter is non-singular). The entropy of the 
output of the channel, i.e., the received signal, will be denoted by //(y). In the noiseless case H(y)=H (x). 
The joint entropy of input and output will be H(xy). Finally there are two conditional entropies H x (y) and 
H y (x), the entropy of the output when the input is known and conversely. Among these quantities we have 
the relations 

H(x, y) = H(x) + H x (y) = H(y) + H y (x). 
All of these entropies can be measured on a per-second or a per-symbol basis. 
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12. Equivocation and Channel Capacity 

If the channel is noisy it is not in general possible to reconstruct the original message o. -the 
signal with certainty by any operation on the received signal E. There are, however, ways of transmuting 
information whfch Ire optima, in combating noise. This is the problem wh.ch ^ 
Suppose there are two possible symbols 0 and 1, and we are transmitting at a rate of ^symbols per 
second'with probabilities po = P, = {■ Thus our source is producing informal,™ a the ra ^o 10M bus 
per second. During transmission the noise introduces errors so that, on the average 1 m 100 « ^recewed 
fnconectlv (a 0 as 1 or 1 as 0). What is the rate of transmission of information? Certamly less than 1000 
oTp" second n e about 1ft of the received symbols are incorrect. Our first impulse rmght be to say 
me ra* So b Us per second, merely subtracting the expected number of errors. Tlus ,s not sadsfect^ 
since it fails to take into account the recipient's lack of knowledge of where the errors occur We may ca^y 
to an extreme case and suppose the noise so great that the received symbols are entirely ^ependen of 
the transmitted symbols. The probability of receiving 1 is \ whatever ^ nM ^^^ m 
Then about half of the received symbols are correct due to chance alone, and we would be giving ^ system 
St remitting 500 bits per second while actually no information is being transmuted at all. Equally 
"good"lnsmission would be obtained by dispensing with the channel entirely and fhppmg a com at the 

^tviSyAe proper correction to apply to the amount of information transmitted is the amount of this 
info™atnwhich P i smissing in the received signal, or f^^^^^^^i 
a sienal of what was actually sent. From our previous discussion of entropy as a measure of uncertainty it 
seems relTab e to use the conditional entropy of the message, knowing the received signal, as a measure 
of to S nformation. This is indeed the proper definition, as we shall see later. Following ttus idea 
L rate oS TSZ transmission, R, would be obtained by subtracting from the rate of products (i.e., the 
entropy of the source) the average rate of conditional entropy. 

R = H(x)-H y (x) 

The conditional entropy H y (x) will, for convenience, be called the equivocation. It measures the average 

is .99, and that a 1 was transmitted is .01. These figures are reversed if a 1 is received. Hence 

H y (x) = -[.991og.99 + 0.01 logO.01] 
= .081 bits/symbol 

or 81 bits per second We may say that the system is transmitting at a rate 1000 - 81 = 919 bits per second 
Z I extrSeTase where a o'is equally likely to be received as a 0 or 1 and similarly for 1, the a postenon 
probabilities are \, \ and 

H y (x) = -[\\og\ + \\og\] 
= 1 bit per symbol 

or 1000 bits per second. The rate of transmission is then 0 as it should be. ,„ ; „ ct;fu 
The following theorem gives a direct intuitive interpretation ofthe equivocation and also serves toju ufy 
it as Ae unfque appropriate measure. We consider a communication system and an observer (or auxihary 
2ele) who can see both what is sent and what is recovered (with errors due to noise), ^observe, - notes 
She errors in the recovered message and transmits data to the receiving point over a "correction channel to 
enable the receiver to correct the errors. The situation is indicated schematically in Fig. 8. 

Theorem 10- If the correction channel has a capacity equal to H y (x) it is possible to so encode the 
correcS^ 

This is not possible if the channel capacity is less than H y (x). 
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CORRECTION DATA 
> — 



OBSERVER 




M 



SOURCE 



TRANSMITTER 



RECEIVER 



CORRECTING 
DEVICE 



Fig. 8 — Schematic diagram of a correction system. 

Roughly then, H y (x) is the amount of additional information that must be supplied per second at the 
receiving point to correct the rece ^, message ; pn ^ , message M < and corresponding original 

To prove the Pf ™^ have produced each 



Hy{x,Z) > Hy{ X ). 



The left-hand side can be expanded to give 



Hy(z) + H yi {x)>H y (x) 
Hyz (x) > H y (x) - Hy{z) > Hy(x) -H(z). 

.fvreidenbfyxasreeo^foftbeso^ 

channel, .hen Ihe right-hand side ts the eqmvocadon lea J^TT^ m0 a „ d 

"g^Kftisis greater .han zero the freqrtencyof erroractmnot be arbtrranly small. 
Example: 

This requires a channel of capacity 

-[plogp + ^logq] 
which is the equivocation of the original system. 
The rate of transmission R can be written in two other forms due to the identities noted above. We have 

R = H(x)-H y (x) 

= H(y)-H x (y) 

= H(x) + H(y)-H(x,y). 
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The first defining expression has already been interpreted as the amount of information sent less the uncer- 
tainty of what was sent. The second measures the amount received less the part of this which is due to noise. 
The third is the sum of the two amounts less the joint entropy and therefore in a sense is the number of bits 
per second common to the two. Thus all three expressions have a certain intuitive significance. 

The capacity C of a noisy channel should be the maximum possible rate of transmission, i.e., the rate 
when the source is properly matched to the channel. We therefore define the channel capacity by 

C = Max(tf(x)-H y (x)) 

where the maximum is with respect to all possible information sources used as input to the channel. If the 
channel is noiseless, H y (x) = 0. The definition is then equivalent to that already given for a noiseless channel 
since the maximum entropy for the channel is its capacity. 

13. The Fundamental Theorem for a Discrete Channel with noise 

It may seem surprising that we should define a definite capacity C for a noisy channel since we can never 
send certain information in such a case. It is clear, however, that by sending the information in a redundant 
form the probability of errors can be reduced. For example, by repeating the message many times and by a 
statistical study of the different received versions of the message the probability of errors could be made very 
small One would expect, however, that to make this probability of errors approach zero, the redundancy 
of the encoding must increase indefinitely, and the rate of transmission therefore approach zero. This* is by 
no means true. If it were, there would not be a very well defined capacity, but only a capacity for a given 
frequency of errors, or a given equivocation; the capacity going down as the error requirements are made 
more stringent. Actually the capacity C defined above has a very definite significance. It is possible to send 
information at the rate C through the channel with as small a frequency of errors or equivocation as desired 
by proper encoding. This statement is not true for any rate greater than C. If an attempt is made to transmit 
at a higher rate than C, say C + R\ , then there will necessarily be an equivocation equal to or greater than the 
excess Ry . Nature takes payment by requiring just that much uncertainty, so that we are not actually getting 
any more than C through correctly. 

The situation is indicated in Fig. 9. The rate of information into the channel is plotted horizontally and 
the equivocation vertically. Any point above the heavy line in the shaded region can be attained and those 
below cannot. The points on the line cannot in general be attained, but there will usually be two points on 
the line that can. 

These results are the main justification for the definition of C and will now be proved. 

Theorem 1 1 : Let a discrete channel have the capacity C and a discrete source the entropy per second H . 
IfH < C there exists a coding system such that the output of the source can be transmitted over the channel 
with an arbitrarily small frequency of errors (or an arbitrarily small equivocation). IfH > C it is possible 
to encode the source so that the equivocation is less thanH-C+e where e is arbitrarily small. There is no 
method of encoding which gives an equivocation less than H-C. 

The method of proving the first part of this theorem is not by exhibiting a coding method having the 
desired properties, but by showing that such a code must exist in a certain group of codes. In fact we will 



Hy( X ) 




C H(x) 
Fig. 9 — The equivocation possible for a given input entropy to a channel. 
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a chnw that this average can be made less than e. If the 

will establish the desired result. 

The capacity C of a noisy channel has been defined as 

C = Max(tf(x)-HyM) 

^ 1 .^^-^^*-----''------**- l * ,ta,,,d "*" 

ud rhe remaining sequent of smUUM.II . M ot lb „„, 2™« members and a low 

2 Similarly the reoerved sequences harm a nrgn proo , 
''rS '« ^ prodoaed o, abour 2™rH> iopobs. The probabilir, of.il olhe, 

* S -5! S525 - «-» -* " abo1 "" io " K, ° " - 

jsss issri: W-s^ — - -* °< caMM for a ,yp,cal 

output. 




2 H{x)T -~ 2 W ^ T 

HIGH PROBABILITY • HIGH PROBABILITY 

MESSAGES ^ — ^ • RECEIVED SIGNALS 

REASONABLE CAUSES 

FOR EACH E • 




REASONABLE EFFECTS # 
FOR EACH M 



Flg . ,0-Schematicrepresentaticnof .he relations between inputs and outputs in a channe.. 

Now suppose we have another source producing ^^^^t^^^V^ 
Z^EXttZSl We will set up ,1s associate in all 
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possible ways (using, however, only the ^^figgX £Tl2 

and average the frequency of errors for this large class « f pos*Me codmg sy ^ 

calculating the frequency of errors for a random asso c auo one message in the set 

T. Suppose a particular output y, is observed. Wham *^ r °3^ rJlW ints . ^ probability of 

of possible causes of y,? There are 2" messages distributed at random P 

a particular point being a message is thus T{R . H (A) 

Theprobabilitythatnoneof thepoints in 

p=[l _ 2 r(«-»W)] 2THyW . 

Now K < H(x) - H y (x) so R - H (*) = -H, W - with „ positive. Consequently 
approaches (as T -» °°) \-2~ Tri . 

Hence probabHi,, of . error approaches aero and feft. ^^^Sc* P e, second 

The second pan of .he theorem ,s easdy shown b, po Png fta « ^ AI „, e 

from the somen, -^-J^SSS «SSidn,» add , This Ital. can also . 

fi = H W = C + ».nd HW -H,W=C + . 

a,,* « positive. This comrade the defini.io. <*J£^™m of « W - »,W • ^ 

14. Discussion 

ThedemonstradonofTh^ 

proofs. An attempt to ffj'^S^T^jS s andcertain limiting situations, no explicit 
generally impracUcal. In fact, apart fron » m '™* ™ found p^y this is no accident but is 
description of a series of approximation to the dart ha been to ' dom nce . 

related to the difficulty of giving an explicit co^^JJ . fa a reasonable way 

A „ approximation to the ^^^^^2-- will not in genera, bring it 
by the no.se, the original can still * Thig h accompli shed at the cost of a certain amount of 

closer to another reasonable signal than the ongmaL 1 m ; is P { ^ ^ 

redundancy in the coding. The ^"^^^T^^^Z^ help if it is utilized at the 
noise structure involved. However any redundancy « «» ^-LcyVno attempt is made to 
receiving point. In particular, if ™ ise For eXamplC> ta 3 

eliminate it in matching to the channe 1 this r ed und ncy «'P f ^ es . This is not done 

telegraph channel one could save about 50% n tune by proper encoding o & * howevefi 

andmo'st of the redundancy of Eng.is h remains in ^a^ymbol, incorrectly 
of allowing considerable noise in the channel. A !^^^ t ^ BUlSm to the ideal in many 

notToo far (in the sense required for the theorem) from a random selection. 
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As in the noiseless case a delay is generally required to approach the ideal encoding. It now has the 
additional function of allowing a large sample of noise to affect the signal before any judgment is made 
at the receiving point as to the original message. Increasing the sample size always sharpens the possible 
statistical assertions. 

The content of Theorem 1 1 and its proof can be formulated in a somewhat different way which exhibits 
the connection with the noiseless case more clearly. Consider the possible signals of duration T and suppose 
a subset of them is selected to be used. Let those in the subset all be used with equal probability, and suppose 
the receiver is constructed to select, as the original signal, the most probable cause from the subset, when a 
perturbed signal is received. We define N(T, q) to be the maximum number of signals we can choose for the 
subset such that the probability of an incorrect interpretation is less than or equal to q. 

Theorem 12: Urn Iglffil^ = C, where C is the channel capacity, provided that q does not equal 0 or 
r->~ T 

1. 

In other words, no matter how we set out limits of reliability, we can distinguish reliably in time T 
enough messages to correspond to about CT bits, when T is sufficiently large. Theorem 12 can be compared 
with the definition of the capacity of a noiseless channel given in Section 1. 

15. Example of a Discrete Channel and its Capacity 

A simple example of a discrete channel is indicated in Fig. 11. There are three possible symbols. The first is 
never affected by noise. The second and third each have probability p of coming through undisturbed, and 
q of being changed into the other of the pair. We have (letting a = -[plogp + qlogq] and P and Q be the 



P 




RECEIVED 
SYMBOLS 



P 

Fig. 1 1 — Example of a discrete channel. 

probabilities of using the first and second symbols) 

H(x) = -P\ogP-2QlogQ 
H y {x) = 2Qa. 

We wish to choose P and Q in such a way as to maximize H(x) - H y (x), subject to the constraint P+2Q = 
Hence we consider 

U = -PlogP-2Q\ogQ-2Qa + \(P+2Q) 

1- logP+A = 0 

2- 21ogG-2a + 2A = 0. 

logP = logQ + a 
P = Qe a = Q(3 



W 
dP 
dU 

Eliminating A 
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The channel capacity is then 0 + 2 

C = log— • 

i o«h n - 1 In the first, 8 = 1 and C = log3, 
H.^ttWto^^^-'"^ ^s. If P = 4. P = 2 and 
which is correct since the channel is then -f^^S a all and acf together like one 
C = log2. Here the second and third ^Xmd and third together with probability 

symbol. Tlte first symbol is used wtth probabthty P - . , and me ^ 
^Thismaybedistnbutedbetweenthem^ The distinction 

2 For intermediate values of p the channel W™VJ™ h as in the nois eless case, 

between the second and third symbols <W«££X * because of its freed om from noise. 
The first symbol is used somewhat more frequently than the 

16 The Channel Capacity in Certain Special Cases 

rate is then given by the maximum of 

- YPiPij^PiPij + 2 p M l0 8P'V 

U ' U 

_ t tn y P . _ ! This leads by the method of Lagrange to the equations, 
where we vary the P, subject to £ n - 1 ■ 1 ms 1Cd » 

v. log Psi -=/i 5 = 1,2...... 

/- t ». tt, P inverse of (if it exists) be hs, so that 
Multiplying by P s and summing on , shows that „ = C Let the mverse ft, 

s.*.w = Then: v - io g 5>> = cis. 



Hence: yp iPir = ex P [-c2^ +2 W°«™] 

° r ' Pi = Vhu expf-C^* + log • 

ofe«ho»>P»''^ w -*"r ; '?r , ^SSSSrfcS symbols, «nd* g iv»b,-X ) > 1 log l >i 

Max[tf(y)-ft(y)] = Maxtf(y) + 5>logP,- 

C = logm + XpilogPi- 
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1/2 

b 

Fig. 12-Examples of discrete channels with the same transition probabilities for each input and for each output. 
In Fig. 12a it would be 

C = log4 - log2 = log2. 
This could be achieved by using only the 1 st and 3d symbols. In Fig. 12b 

C = log4-2log3-^log6 
= lbg4 - log3 - \ log2 
= logi23. 



In Fig. 12c we have 



C = log3 - \ log2 - \ log3 - \ log6 
3 

= to8 .i.i 



2Z3365 

Suppose the symbols fall into several groups such that the noise never causes a symbol in one group to 
be mistaken for a symbol in another group. Let the capacity for the nth group be C„ (in bits per second) 
when we use only the symbols in this group. Then it is easily shown that, for best use of the entire set, the 
total probability P„ of all symbols in the nth group should be 



2 Cn 



Within a group the probability is distributed just as it would be if these were the only symbols being used. 
The channel capacity is 

C = log£2 c \ 
17. An Example of Efficient Coding 

The following example, although somewhat unrealistic, is a case in which exact matching to a noisy channel 
^Sl There are two channel symbols, 0 and 1, and the noise affects them in blocks of seven symbols. 
A Sto* of seven is either transmitted without error, or exactly one symbol of the seven is incorrect. These 
eight possibilities are equally likely. We have 

C = M*x[H(y)-H x (y)] 

= 7[7 + |log£] 
= | bits/symbol. 

An efficient code, allowing complete correction of errors and transmitting at the rate C, is the following 
(found by a method due to R. Hamming): 
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. y y, v. of these X 3 , X 5 , X 6 and X 7 are message symbols and 
chosen aAitiajily tythe Murce/The other diree are reJundant sad calculated as follows: 

X 4 is chosen to make a = X 4 + X 5 + X 6 +X 7 even 
X2 - " " " ^ = X2 + X 3 + X 6 + X 7 
Y j " " " " 7 = Xj + X3 + X5 + X7 

• a ft a „H^ fl rf calculated and if even called zero, if odd called one. The 
APPENDIX 1 

THE GROWTH OF THE NUMBER OF BLOCKS OF SYMBOLS WITH A FINITE STATE CONDITION 
Let Ni{L) be the number of blocks of symbols of length L ending in state , Then we have 

where *»,*„...,«! are the length of the symbols which may be chosen in state i and lead to state These 
S d^V» and the behavior as L -> ~ must be of the type 



Substituting in the difference equation 



Nj=AjW 



or 



X(lW-^ , -«5, 7 )A i = 0. 

For this to be possible the determinant 

d(w) = M = |S^ ;) -^| 

must vanish and this determines W, which is, of course, the largest real root of D = 0. 
The quantity C is then given by 

and we also note that the same growth proves result if we require that all blocks start in the same (arbi- 
trarily chosen) state. 

APPENDIX 2 
Derivation ofH = -Sp/logp/ 
LetHfi i I) = A(n). From condition (3) we can decompose a choice from y» equally likely possi- 
bilities into a series of m choices from s equally likely possibilities and obtain 

A{s m )=mA(s). 
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Similar,y A(f n )=nA(r). 
We can choose n arbitrarily large and find an m to satisfy 

f <t" <s (m+,) . 

Thus, taking logarithms and dividing by nlogs, 

n " logs" n n l« lo 8 sl 
where e is arbitrarily small. Now from the monotonic property of A(n), 

A(^)<A{t n )<A(^ +x ) 
mA(s)<nA(t)<{m+l)A(s). 

Hence, dividing by nA(s), ^ |m_M|< e 

|M0_!2iH<2e ,*(,) = Klogf 
U(s) logs I 

J^logX"; = H(Pi>- -,P n ) + ^' l0 S ni - 



Hence 



^A^Xp/logSnz-XPilog"'-] 



= -^Sp/log^ r = -^SPi 1 °gP'- 

of convenience and amounts to the choice of a unit of measure. 

APPENDIX 3 
Theorems on Ergodic Sources 

. with p ^ 0 to anv other along a path of probability p > 0, the system is 
If it is possible to go from any state with P > 0 m an otoera g P ^ ^ ^ ^ ^ ^ 

ergodic and the strong law of large numbers «teijJ0 , ^ probab ility of being at 

the network is traversed in a long '^/^^^S^^ of percentage error ±5 in 

{P iPij ±6)N. 

Hence nearly all sequences have a probability p given by 

n(P iPi j±S)N 
Pij 
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and !2i£ is limited by 

!^ = £(P iW ±5)logp,y 

possible range of values of p in Theorem 3. 
In the mixed (not ergodic) case if 

and the entropies of the components are H, > «2 > ■ ■ ■ > ft we have the 
Theorem: Lim ^ = <?( 9 ) is a decreasing step function. 



5-1 



<p{q) = ft in the intervaJ X ai < q K £ a '' 



F N = NG N -(N-l)G N -\ 



, ,. x, • /- 1 V P Hence Gw > Fv and Gat monotonic decreasing. Also they 
and summing this for all N gives G N = jj 2A Hence ° w - ™ _ 
must approach the same limit. By using Theorem 3 we see that Urn G N - H. 

APPENDIX 4 

MAXIMIZING THE RATE FOR A SYSTEM OF CONSTRAINTS 



Pi)' 

The rate in question is . . 



Let I ■ = S # Evidently for a maximum p» = *exp$ . The constraints on maximization are Ifl = 

i. S;-py = i. i flU - *y) = °- Hence we maxinuze 

?ui_ _ MP.(i +iogPu)+NPitij | ^ + f ,, ^.^.p. = n . 
3/»y M " 
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Solving for py 

Pij^AiBjD-t*. 

Since 

j J 
BjD- li > 

The correct value of D is the capacity C and the Bj are solutions of 

Bi^BjC^ 



for then 



or 



So that if A; satisfy 



^Bi Bj' 



Pi = Biji. 



Both the sets of equations for B, and 7, can be satisfied since C is such that 

|C-<<;_<Sij|=0. 

In this case the rate is « . » 

but „ 

^PiPijOogBj - logB,) = 2P;logB; - XfilogB, = 0 

j 

Hence the rate is C and as this could never be exceeded this is the maximum, justifying the assumed solution. 
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PART HI: MATHEMATICAL PRELIMINARIES 

In mis final M of the paper we consider Jj-jj — S 
continuously variable, in contrast with the discrete ^^^^case by dividing the continuum 
continuouscasecanbeobtainedm^ 
ofmessages and signals ^alargebutfinnenu^ 

involved on a discrete basis. As the size of * e J e ^^^~ er , few new effects that appear and also 
limits the proper values for the continuous case^ There are . however a t ^ ^ 

a general change of emphasis in the direction of ^^^^ gre atesi generality, or with 
We will not attempt, in the continuous case to obtain our r« ul ts wr«i g 6 ^ 

can be justified in all cases of practical interest. 

18. Sets and Ensembles of Functions 
We shall have to deal in the continuous case with 

functions, as the name implies, is merely a in the set, or implicitly by 

1. Thesetof functions: ^ =sin(f + e) . 

Each particular value of 9 determines a particular function in the set. 

2. The set of all functions of time containing no frequencies over W cycles per second. 

3. The set of all functions limited in band to W and in amplitude to A. 
4 The set of all English speech signals as functions of time. 

^^^^^^^^ 

/„(r) = sin(r + 0), 

we may give a probability distribution for 9, P(9) . The set then becomes an ensemble. 
Some further examples of ensembles of functions are: 

1. A finite set of functions f k (t) (*= 1,2,..., ») with the probability of /.being p k . 

2. A finite dimensional family of functions 

f(a\,at,...,ct n \i) 

with a probability distribution on the parameters a,-: 

p(a\,...,a n ). 

For example we could consider the ensemble defined by 

71 

/(«,,. ..,«.,«! ».:<) = g«nM«+*> 

(from 0 to 2tt) and independently. 
■ In nuthemaoca. termino.ogy ft. functions be.ong lo a measure space whose to*, measure is unity. 



32 



3. The ensemble +- s in7 r(2Wf - n) 

/(«..') = X a " *(2Wt-n) 

i /('+<•) 

tapulse o, sho, noise where all me impute are .(renocal. 

in ordinary use. , 
any fixed amount in time. The ensemble 

/,(*) = sinO+0) 

=«&»(' 

= sin(/ + </?) 

* , * ntn7ir Each function has changed but the ensemble as a whole is 
with f distributed uniformly from ^^^L. are also stationary. 

inv t — - « ?r of the funct,ons in the set a 

probability different from 0 and 1 which is stationary. The ensemble 

sin(f + 0) 

lations. On the other hand the ensemble 

asin\f-rtv 

a between 0 and 1 for example is d 5 haps be considered so. If an ensemble is 

Of the examples given, 3 and 4 are ^ergodic, _ and 3 may PP We More predsely It „ 

of any of the functions in the set 

^rep^fioncanheuseaasa-efi^ 
Umiting operations than do definitions that have been ^» J^^TEy confinuous spectrum as contrasted wdh a porn, 
literature, is perhaps somewhat unfortunate. In ophcs wtatt ^ B , s £ ctro m flat with frequency). 
™Lmim or r S pectrum which is flat with wavelength (which .s aot dean : h ' somewhat different formulations 

SP ^'is tne toous ergodic theorem or rather one aspect fJ^^^Z^^cz and others. Tire Uterature on 
by BMrff von Neumann, and Koopman, and subsequenUy genend^d J ise ^ eral formu latious ; e.g.. 

ereSc theory is quite extensive and the reader is referred to the > papefi of Statistics ^ Probability." Jtanrf 

of Mathematics and Physics, v. XIII, No. 1, 1934, N. W.ener, h» 
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Just as we may perform various operations on numbers or functions to obtain new numbers or funcUons 
we can perform operations on ensembles to obtain new ensembles. Suppose, or ex we have an 
Tsemble of funds f a (t) and an operator T which gives for each funcUon f a (t) a resultmg funcUon 

8a(t): ga (t) = Tf a (t). 

Probability measureisdefinedfor the set^^ 

^bset of L g a (t) functions is equal to that of the subset of the f a (t) functions wh.ch produce members of 
AeSven subtt of g functions ler the operation T. Physically this corresponds to passing *e ensemWe 
through some device, for example, a filter, a rectifier or a modulator. The output funcUons of the dev.ce 

^A^Zofo^T will be called invariant if shifting the input merely shifts the output, i.e., if 

g a (t) = Tf a (t) 

implies , „ , . . 

g a (t+h) = Tf a (t+h) 

for all f (t) and all r,. It is easily shown (see Appendix 5 that if T is invariant and the input ensemble is 
Sti^ *en the output ensemble is stationary. Likewise if the input is ergodic the output w,l. also be 

^Afiher or a rectifier is invariant under all time translations. The operation of modulation is not since Mite 
carrier phase gives a certain time structure. However, modulation is invariant under all translates which 

^ tSESSSSZZSL Nation between the invariance of physical devices under time ' 
tranlfions andFourier theory.' He has shown, in fact, that if a device is linear as well as mvanant Founer 
analvsis is then the appropriate mathematical tool for dealing with the problem. 

An ens mble of funds is the appropriate mathematical representation of the mess ag« ^ 
a continuous source (for example, speech), of the signals produced by a transmitter, and of the perturbing 

onT^ 

notfor a particular speech function and still less for a sine wave, but for the ensemble of speech functions. 

19. Band Limited Ensembles of Functions 

If a function of time f{t) is limited to the band from 0 to W cycles per second it is com P lete ^ e ™ n ^ 
by giving its ordinate* at a series of discrete points spaced £ seconds apart m the manner indicated by the 

following result. 5 

Theorem 13: Let /(f) contain no frequencies over W. Then 

f(t)-Yx sin * {2Wt - n) 



where 



Wnunication theory is heavUy indebted le aner for much of Us ^^ n ^ fiK ^^ f S 

unear preu^uuu 5 * g ^ ith me neml pro blems of communication and control. 

*"* S^tTS^ see JLuJs paper "Co—on in «he Presence of Noise" pubhshed ,n 

the Proceedings of the Institute of Radio Engineers, v. 37, No. 1, Jan., 1949, pp. 10-21. 
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function corresponds to precisely .one .point and each toon ^x, outride ihi. 

a band W and duration T correspond to points in J.^"^^ in mis space . For example, 

P (, 1) ... > x n )inthecone S ponding«d—o 

the 2HV coordinates m a given interval 7 to ««cL of the ensemble for intervals of 

and the probability distribution p(*i , . . • ,*«) to give the siausnwi 
that duration. 

20. ENTROPY OF A CONTINUOUS DISTRIBUTION 
The entropy of a discrete set of probabilities p, , . . . , Pn has been defined as: 

f/ = -£p/10gPi- 

In an analogous manner we define the entropy of a continuous distribution with the density distribution 

function p(x) by: „ 

H = - I p{x)\ogp{x)dx. 

With an n dimensional distribution p(*i , . . • ,x„) we have 

fj = - J ■ ■ J p(x\ ,x n ) logp(*i , • • • ,x n )dx\ ■ ■ ■ dx„. 

If W e have two arguments , and y (which may themselves be multidimensional) the joint and conditional 
entropies of p{x 7 y) are given by 

■H(xj) = -jf p{x,y)\ogp{x 7 y)dxdy 

and 

H y (x) = -JJp(x,y)\og^-dxdy 

where 

p{x) = f p{x,y)dy 
p(y) = J P(x,y)dx. 

The entropies of continuous distributions have most (but not a.l) of the properties of the discrete case. 
In particular we have the following: 

1. If x is limited to a certain volume v in its space, then H(x) is a maximum and *mal to logv when ,(,) 
is constant (1/v) in the volume. 
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2. With any two variables x, y we have 

H(x,y)<H(x) + H(y) 

.» it <M on,y if) , .d , - ■ ' »M - <*« *-* «" ' 

set of points of probability zero). 

3. Consider a generalized averaging operation of the following type: 

p'(y) = J a{x,y)p{x)dx 

^ Ja( X ,y)dx = Ja { x,y)dy=l, a{x,y)>0. 

Then the entropy of the averaged distribution , « is equal to or greater than mat of the original 
distribution p(x). 
4. We have 

H(x;y) = H W + ftGO = + 

and 

ft(y)<ff(y)- 

H(x) = - J p(x)logp(x)dx 

With a ^Jp{x)x 1 dx and l = fp(x)dx 

as constraints. This requires, by the calculus of variations, maximizing 

J [-p(x) \ogp(x) + Xp (x)x 2 + «K*)] dx - 

The condition for this is _ t _ logp(jc) + A* 2 + p = 0 

and consequently (adjusting the constants to satisfy the constraints) 

Similarly in n dimensions, suppose the second order moments of pfr *„) are fixed at A* 

Aij = /•••/ **;p(*i • • ■ ■ dxi " dXn - 

Thenmemaximumentropyoccurs(byasimilarcalcu,ation)w^ 
Gaussian distribution with the second order moments Ay. 
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6. The entropy of a one-dimensional Gaussian distribution whose standard deviation is a is given by 

tf(x) = log\/27r«r. 

This is calculated as follows: 

JC 2 

-logp(x) = log\/27r<7+ ^2 
H(jc) = - f p{x)\ogp{x)dx 

= y p(x)tog>/2ir<Tdx + j pW^2^ 

= logV27rcr+^2 
= log\/27rcr + logv^ 

= logV2?T€(7. 

Similarly the n dimensional Gaussian distribution with associated quadratic form a u is given by 

p(« = 

and the entropy can be calculated as 

H = log(27r^ 2 |a j7 r^ 

where |a;j| is the determinant whose elements are ay. 
7. If x is limited to a half line (p(x) = 0 for x < 0) and the first moment of x is fixed at a: 

JO 

then the maximum entropy occurs when 
and is equal to logea. 

8 There is one important difference between the continuous and discrete entropies. In the discrete case 
JS^ZSTin an absolute way the randomness of the chance variable. In the contmuou 
^ S«nt is relative to the coordinate system. If we change coordinates the entropy w.ll 
I T^TchTge. In fact if we change to coordinates * • • -y n the new entropy » gwen by 

W (v) = /-/ P (x ll ... 1 x n y(^)logp(x 1) ... ) x n )7^)dyi-^ 

where 7(f) is the Jacobian of the coordinate transformation. On expanding the logarithm and chang- 
ing the variables to x\ ■■ ■ x„, we obtain: 

H{y) = H(x)-f-Jp{xi,...,x„)\ogj{^)dx l ...dx n . 
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Thus the new entropy is the old entropy less the expected logarithm of the Jacobian. In the continuous 
case the entropy can be considered a measure of randomness relative to an assumed standard, namely 
the coordinate system chosen with each small volume element dx\ • • dx n given equal weight. When 
we change the coordinate system the entropy in the new system measures the randomness when equal 
volume elements dy\ • • • dy n in the new system are given equal weight. 

In spite of this dependence on the coordinate system the entropy concept is as important in the con- 
tinuous case as the discrete case. This is due to the fact that the derived concepts of information rate 
and channel capacity depend on the difference of two entropies and this difference does not depend 
on the coordinate frame, each of the two terms being changed by the same amount. 

The entropy of a continuous distribution can be negative. The scale of measurements sets an arbitrary 
zero corresponding to a uniform distribution over a unit volume. A distribution which is more confined 
than this has less entropy and will be negative. The rates and capacities will, however, always be non- 
negative. 

9. A particular case of changing coordinates is the linear transformation 

yj = Say**- 

In this case the Jacobian is simply the determinant \aij\~ 1 and 

H(y) = H(x) + log\a ij \. 

In the case of a rotation of coordinates (or any measure preserving transformation) 7=1 and H(y) = 
H(x). 

21. Entropy of an Ensemble of Functions 
Consider an ergodic ensemble of functions limited to a certain band of width W cycles per second. Let 

p(x u .,.,x n ) 

be the density distribution function for amplitudes x\ , . . . , x n at n successive sample points. We define the 
entropy of the ensemble per degree of freedom by 

H' = -Um- A ■ f p{x\,..^x n )\ogp(x u ...,x n )dx\ ,,.dx n . 
n->~ n J J 

We may also define an entropy H per second by dividing, not by n, but by the time T in seconds for n 
samples. Since n = 2TW, H = 2WH'. 

With white thermal noise p is Gaussian and we have 

H f = \ogy/2neN 1 
H = W\og2neN. 

For a given average power N, white noise has the maximum possible entropy. This follows from the 
maximizing properties of the Gaussian distribution noted above. 

The entropy for a continuous stochastic process has many properties analogous to that for discrete pro- 
cesses. In the discrete case the entropy was related to the logarithm of the probability of long sequences, 
and to the number of reasonably probable sequences of long length. In the continuous case it is related in 
a similar fashion to the logarithm of the probability density for a long series of samples, and the volume of 
reasonably high probability in the function space. 

More precisely, if we assume p(x\ , . . . continuous in all the for all n, then for sufficiently large n 

\ l ^-H'\« 
I n I 



38 



for all choices of (x\ , . . . ,*„) apart from a set whose total probability is less than £, with 8 and e arbitrarily 
small. This follows form the ergodic property if we divide the space into a large number of small cells. 

The relation of H to volume can be stated as follows: Under the same assumptions consider the n 
dimensional space corresponding to p(xi,...,x„). Let V n (q) be the smallest volume in this space which 
includes in its interior a total probability q. Then 

n 

provided q does not equal 0 or 1. 

These results show that for large n there is a rather well-defined volume (at least in the logarithmic sense) 
of high probability, and that within this volume the probability density is relatively uniform (again in the 
logarithmic sense). 

In the white noise case the distribution function is given by 

^''•••' Xn)= (2^ eXp -2}v^- 

Since this depends only on £jc? the surfaces of equal probability density are spheres and the entire distri- 
bution has spherical symmetry. The region of h igh proba bility is a sphere of radius yfnN. As n -> «> the 
probability of being outside a sphere of radius y/n(N + e) approaches zero and i times the logarithm of the 
volume of the sphere approaches \og\f2ixeN. 

In the continuous case it is convenient to work not with the entropy H of an ensemble but with a derived 
quantity which we will call the entropy power. This is defined as the power in a white noise limited to the 
same band as the original ensemble and having the same entropy. In other words if H* is the entropy of an 
ensemble its entropy power is 

N\ = -i-exp2//'.. 
lire v 

In the geometrical picture this amounts to measuring the high probability volume by the squared radius of a 
sphere having the same volume. Since white noise has the maximum entropy for a given power, the entropy 
power of any noise is less than or equal to its actual power. 

22. Entropy Loss in Linear Filters 

Theorem 14: If an ensemble having an entropy H\ per degree of freedom in band W is passed through a 
filter with characteristic Y(f) the output ensemble has an entropy 



H 2 =H { + ± [ \og\Y(f)\ 2 df. 
W Jw 



The operation of the filter is essentially a linear transformation of coordinates. If we think of the different 
frequency components as the original coordinate system, the new frequency components are merely the old 
ones multiplied by factors. The coordinate transformation matrix is thus essentially diagonalized in terms 
of these coordinates. The Jacobian of the transformation is (for n sine and n cosine components) 

1=1 

where the f\ are equally spaced through the band W. This becomes in the limit 

exp i^ l w ^\nf)\ 2 df. 

Since J is constant its average value is the same quantity and applying the theorem on the change of entropy 
with a change of coordinates, the result follows. We may also phrase it in terms of the entropy power. Thus 
if the entropy power of the first ensemble is N\ that of the second is 



J w \og\Y{f)\ 2 df. 
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TABLE I 



GAIN 


ENTROPY 
POWER 
FACTOR 


ENTROPY 
POWER GAIN 
IN DECIBELS 


IMPULSE RESPONSE 


1 

1 >- 

0 


w 1 


1 

? 


-8.69 


sin 2 (r/2) 
f 2 /2 


1 

j 1-uj 2 --->" 
0 


u 1 


/2\4 

(;) 


-5.33 


^Tsinr cosrl 

2 [— -— J 


1 

0 


w 1 


0.411 


—3.87 


[cosr-l cosr sinrl 
6 [ r\ 2/ 2 + .|3 J 


1 

Vi-w 2 — 

0 


w 1 




— z.o / 


2 t 

V 


i 

0 


W 1 


1 

e 2a 


-8.69q 


-L fcos( 1 - a)r - cos/1 
af 2 1 J 



The final entropy power is the initial entropy power multiplied by the geometric mean gain of the filter. If 
the gain is measured in db, then the output entropy power will be increased by the arithmetic mean db gain 
over W. 

In Table I the entropy power loss has been calculated (and also expressed in db) for a number of ideal 
gain characteristics. The impulsive responses of these filters are also given for W = 27r, with phase assumed 
to be 0. 

The entropy loss for many other cases can be obtained from these results. For example the entropy 
power factor l /e 2 for the first case also applies to any gain characteristic obtain from 1 — a; by a measure 
preserving transformation of the u axis. In particular a linearly increasing gain G{lj) = a;, or a "saw tooth" 
characteristic between 0 and 1 have the same entropy loss. The reciprocal gain has the reciprocal factor. 
Thus 1 /w has the factor e 2 . Raising the gain to any power raises the factor to this power. 

23. Entropy of a Sum of Two Ensembles 

If we have two ensembles of functions f Q (t) and gp(t) we can form a new ensemble by "addition." Suppose 
the first ensemble has the probability density function p(x\ , . . . ,*„) and the second q(x\ , . . . ,*„). Then the 
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density function for the sum is given by the convolution: 

r{x u ...,*„) = J-J p(y\,--->yn)q(x\-y\>.--,x n -y n )dyi-~dy n . 

Physically this corresponds to adding the noises or signals represented by the original ensembles of func- 
tions. 

The following result is derived in Appendix 6. 

Theorem 75: Let the average power of two ensembles be N\ and N 2 and let their entropy powers be N \ 
and N 2 . Then the entropy power of the sum, N 3 , is bounded by 

Ni+N 2 <N3<N\+N 2 . 

White Gaussian noise has the peculiar property that it can absorb any other noise or signal ensemble 
which may be added to it with a resultant entropy power approximately equal to the sum of the white noise 
power and the signal power (measured from the average signal value, which is normally zero), provided the 
signal power is small, in a certain sense, compared to noise. 

Consider the function space associated with these ensembles having n dimensions. The white noise 
corresponds to the spherical Gaussian distribution in this space. The signal ensemble corresponds to another 
probability distribution, not necessarily Gaussian or spherical. Let the second moments of this distribution 
about its center of gravity be a\ y That is, if p(x\ , . . . ,x n ) is the density distribution function 

a/j = / '/ P( x i- a d( x j- a j) dx i"' dx n 

where the a; are the coordinates of the center of gravity. Now ay is a positive definite quadratic form, and 
we can rotate our coordinate system to align it with the principal directions of this form. a\ } is then reduced 
to diagonal form We require that each bn be small compared to N, the squared radius of the spherical 
distribution. 

In this case the convolution of the noise and signal produce approximately a Gaussian distribution whose 
corresponding quadratic form is 

N + b u . 

The entropy power of this distribution is 

[U( N + b 4 Un 

or approximately 

= N + -2bii. 
n ^ 

The last term is the signal power, while the first is the noise power. 

PART IV: THE CONTINUOUS CHANNEL 

24. The Capacity of a Continuous Channel 

In a continuous channel the input or transmitted signals will be continuous functions of time /(f) belonging 
to a certain set, and the output or received signals will be perturbed versions of these. We will consider 
only the case where both transmitted and received signals are limited to a certain band W. They can then 
be specified, for a time 7\ by 2TW numbers, and their statistical structure by finite dimensional distribution 
functions. Thus the statistics of the transmitted signal will be determined by 

P(x u ... )Xn )=P{x) 
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and those of the noise by the conditional probability distribution 

The rate of transmission of information for a continuous channel is defined in a way analogous to that 
for a discrete channel, namely 

R = H{x)-H y {x) 

where H(x) is the entropy of the input and H y (x) the equivocation. The channel capacity C is defined as the 
maximum ofR when we vary the input over all possible ensembles. This means that in a finite dimensional 
approximation we must vary P(x) = P(x\ , . . . ,x„) and maximize 

P(x,y) 



- J P(x)\ogP{x)dx + If P( x ,y)log-^dxdy. 



This can be written 

using the fact that Jj P(x,y) \ogP(x)dxdy = J P(x) \ogP{x)dx. The channel capacity is thus expressed as 
follows: 

It is obvious in this form that R and C are independent of the coordinate system since the numerator 
P(x y) 

and denominator in log '\ ,\ will be multiplied by the same factors when x and y are transformed in 

P(x)P{y) . 
any one-to-one way. This integral expression for C is more general than H (x) - H y (x). Properly interpreted 
(see Appendix 7) it will always exist while H(x) - H y (x) may assume an indeterminate form oo - oo in some 
cases. This occurs, for example, if x is limited to a surface of fewer dimensions than n in its n dimensional 
approximation. 

If the logarithmic base used in computing H(x) and H y (x) is two then C is the maximum number of 
binary digits that can be sent per second over the channel with arbitrarily small equivocation, just as in 
the discrete case. This can be seen physically by dividing the space of signals into a large number of 
small cells, sufficiently small so that the probability density P x (y) of signal x being perturbed to point y is 
substantially constant over a cell (either of x or y). If the cells are considered as distinct points the situation is 
essentially the same as a discrete channel and the proofs used there will apply. But it is clear physically that 
this quantizing of the volume into individual points cannot in any practical situation alter the final answer 
significantly, provided the regions are sufficiently small. Thus the capacity will be the limit of the capacities 
for the discrete subdivisions and this is just the continuous capacity defined above. 

On the mathematical side it can be shown first (see Appendix 7) that if u is the message, x is the signal, 
y is the received signal (perturbed by noise) and v is the recovered message then 

H(x)-H y (x)>H(u)-H v (u) 

regardless of what operations are performed on u to obtain x or on y to obtain v. Thus no matter how we 
encode the binary digits to obtain the signal, or how we decode the received signal to recover the message, 
the discrete rate for the binary digits does not exceed the channel capacity we have defined. On the other 
hand, it is possible under very general conditions to find a coding system for transmitting binary digits at the 
rate C with as small an equivocation or frequency of errors as desired. This is true, for example, if, when we 
take a finite dimensional approximating space for the signal functions, P{x,y) is continuous in both x and y 
except at a set of points of probability zero. 

An important special case occurs when the noise is added to the signal and is independent of it (in the 
probability sense). Then P x (y) is a function only of the difference n = (y -*), 

P x (y) = Q(y-x) 
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and we can assign a definite entropy to the noise (independent of the statistics of the signal), namely the 
entropy of the distribution Q(n). This entropy will be denoted by H (n). 

Theorem 16: If the signal and noise are independent and the received signal is the sum of the transmitted 
signal and the noise then the rate of transmission is 

R = H(y)-H(n), 

i.e., the entropy of the received signal less the entropy of the noise. The channel capacity is 

C = MaxH(y)-H(n). 

P(x) 

We have, since y = x + n: 

H(x,y)=H(x,n). 
Expanding the left side and using the fact that x and n are independent 

H(y) + H y (x)=H(x) + H(n). 

R = H(x) - H y (x) = H(y) - H(n). 

Since H(n) is independent of P(x) y maximizing R requires maximizing tf(y), the entropy of the received 
signal. If there are certain constraints on the ensemble of transmitted signals, the entropy of the received 
signal must be maximized subject to these constraints. 

25. Channel Capacityiwith an Average Power Limitation 

A simple application of Theorem 16 is the case when the noise is a white thermal noise and the transmitted 
signals are limited to a certain average power P. Then the received signals have an average power P + N 
where N is the average noise power. The maximum entropy for the received signals occurs when they also 
form a white noise ensemble since this is the greatest possible entropy for a power P + N and can be obtained 
by a suitable choice of transmitted signals, namely if they form a white, noise ensemble of power P. The 
entropy (per second) of the received ensemble is then 

H(y) = Wlog2ire(P + N) r 

and the noise entropy is 

H{n)=W\og2neN. 

The channel capacity is 

P + N 

C = H(y)-H(n) = W\og — . 
Summarizing we have the following: 

Theorem 17: The capacity of a channel of band W perturbed by white thermal noise power N when the 
average transmitter power is limited to P is given by 

This means that by sufficiently involved encoding systems we can transmit binary digits at the rate 
\y j 0 g 2 bits per second, with arbitrarily small frequency of errors. It is not possible to transmit at a 
higher rate by any encoding system without a definite positive frequency of errors. 

To approximate this limiting rate of transmission the transmitted signals must approximate, in statistical 
properties, a white noise. 6 A system which approaches the ideal rate may be described as follows: Let 

6 This and other properties of the white noise case are discussed from the geometrical point of view in "Communication in the 
Presence of Noise," loc. cit. 
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M = 2 s samples of white noise be constructed each of duration 7. These are assigned binary numbers from 
0 to M - 1. At the transmitter the message sequences are broken up into groups of s and for each group 
the corresponding noise sample is transmitted as the signal. At the receiver the M samples are known and 
the actual received signal (perturbed by noise) is compared with each of them. The sample which has the 
least R.M.S. discrepancy from the received signal is chosen as the transmitted signal and the corresponding 
binary number reconstructed. This process amounts to choosing the most probable (a posteriori) signal. 
The number M of noise samples used will depend on the tolerable frequency e of errors, but for almost all 
selections of samples we have 

. T . logAf (c,T) ■ P + N 
Lim Lim & ^ = W log — — , 

so that no matter how small e is chosen, we can, by taking T sufficiently large, transmit as near as we wish 
P + N 

to TW log — — binary digits in the time T. 

Formulas similar to C = Wlog^p for the white noise case have been developed independently 
by several other writers, although with somewhat different interpretations. We may mention the work of 
N. Wiener, 7 W. G. Tuller, 8 and H. Sullivan in this connection. 

In the case of an arbitrary perturbing noise (not necessarily white thermal noise) it does not appear that 
the maximizing problem involved in determining the channel capacity C can be solved explicitly. However, 
upper and lower bounds can be set for C in terms of the average noise power N the noise entropy power N\. 
These bounds are sufficiently close together in most practical cases to furnish a satisfactory solution to the 
problem. 

Theorem 18: The capacity of a channel of band W perturbed by an arbitrary noise is bounded by the 
inequalities . . mi 

wlog ^<C<W.og^ 

where 

P = average transmitter power 
N == average noise power 
N\ = entropy power of the noise. 

Here again the average power of the perturbed signals will be P + N. The maximum entropy for this 
power would occur if the received signal were white noise and would be W\og2ne(P + N). It may not 
be possible to achieve this; i.e., there may not be any ensemble of transmitted signals which, added to the 
perturbing noise, produce a white thermal noise at the receiver, but at least this sets an upper bound to H(y). 
We have, therefore 

C = MaxH(y)-H(n) 

< W log 2-ne{P + N) - W log 2neN\ . 

This is the upper limit given in the theorem. The lower limit can be obtained by considering the rate if we 
make the transmitted signal a white noise, of power P. In this case the entropy power of the received signal 
must be at least as great as that of a white noise of power P + N\ since we have shown in in a previous 
theorem that the entropy power of the sum of two ensembles is greater than or equal to the sum of the 
individual entropy powers. Hence 

Maxtf(y) > W\og27re(P + Ni) 

7 Cybernetics, loc. tit. . 
8t Theoretical Limitations on the Rate of Transmission of Information," Proceedings of the Institute of Radio Engineers, v. J/, 

No. 5, May, 1949, pp. 468-78. 
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and 



C > W \og2ire(P + N\)-W log27reJVi 
As P increases, the upper and lower bounds approach each other, so we have as an asymptotic rate 

If the noise is itself white, N = N\ and the result reduces to the formula proved previously: 



C = Wlog(l + £). 



If the noise is Gaussian but with a spectrum which is not necessarily flat, N\ is the geometric mean of 
the noise power over the various frequencies in the band W . Thus 



W Jw 



where N(f) is the noise power at frequency /. 

Theorem 19: If we set the capacity for a given transmitter power P equal to 

N\ 

then 77 is monotonic decreasing as P increases and approaches Oasa limit. 
Suppose that for a given power P\ the channel capacity is 

Pi+N-m 



Wlog- 



This means that the best signal distribution, say p(x) y when added to the noise distribution q(x), gives a 
received distribution r{y) whose entropy power is (P\ + N — t)\). Let us increase the power to P\ + AP by 
adding a white noise of power AP to the signal. The entropy of the received signal is now at least 

H(y)=W\og2Tre(Pi+N-r}i+AP) 

by application of the theorem on the minimum entropy power of a sum. Hence, since we can attain the 
H indicated, the entropy of the maximizing distribution must be at least as great and 77 must be monotonic 
decreasing. To show that 77 -> 0 as P -» <» consider a signal which is white noise with a large P. Whatever 
the perturbing noise, the received signal will be approximately a white noise, if P is sufficiently large, in the 
sense of having an entropy power approaching P + N. 

26. The Channel Capacity with a Peak Power Limitation 

In some applications the transmitter is limited not by the average power output but by the peak instantaneous 
power. The problem of calculating the channel capacity is then that of maximizing (by variation of the 
ensemble of transmitted symbols) 

H(y)-H(n) 

subject to the constraint that all the functions f(t) in the ensemble be less than or equal to VS, say, for all 
r. A constraint of this type does not work out as well mathematically as the average power limitation. The 

most we have obtained for this case is a lower bound valid for all an "asymptotic" upper bound (valid 

N 

S S 
for large — ) and an asymptotic value of C for — small. 
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Theorem 20: The channel capacity C for a band W perturbed by white thermal noise of power N is 
bounded by 

S 

where S is the peak allowed transmitter power. For sufficiently large — 

JLc 4. pj 
C<W\og " e N (1+c) 

5 

where e is arbitrarily small. As — -> 0 (and provided the band W starts at 0) 



:/wiog(i + £) 



1. 



We wish to maximize the entropy of the received signal. If — is large this will occur very nearly when 

we maximize the entropy of the transmitted ensemble. 

The asymptotic upper bound is obtained by relaxing the conditions on the ensemble. Let us suppose that 
the power is limited to S not at every instant of time, but only at the sample points. The maximum entropy of 
the transmitted ensemble under these weakened conditions is certainly greater than or equal to that under the 
original conditions. This altered problem can be solved easily. The maximum entropy occurs if the different 
samples are independent and have a distribution function which is constant from — \/S to +y/S. The entropy 
can be calculated as 

Wlog4S. 

The received signal will then have an entropy less than 

Wlog(4S + 27r<>A0(l+e) 

with e -¥ 0 as -» oo and the channel capacity is obtained by subtracting the entropy of the white noise, 

W\og2neN: 

—S + N 

W log(45 + 2neN) ( 1+ e) - W \og(2ireN) = W log ss— — ( 1 + e) . 

This is the desired upper bound to the channel capacity. 

To obtain a lower bound consider the same ensemble of functions. Let these functions be passed through 
an ideal filter with a triangular transfer characteristic. The gain is to be unity at frequency 0 and decline 
linearly down to gain 0 at frequency W. We first show that the output functions of the filter have a peak 

. x ^. . sin27rWr 

power limitation S at all times (not just the sample points). First we note that a pulse _ . going into 

the filter produces 

1 sin 2 7rWr 

2 (ttW/) 2 

in the output. This function is never negative. The input function (in the general case) can be thought of as 
the sum of a series of shifted functions 

sin27r\V/ 
a 2irWt 

where a> the amplitude of the sample, is not greater than y/S. Hence the output is the sum of shifted functions 
of the non-negative form above with the same coefficients. These functions being non-negative, the greatest 
positive value for any t is obtained when all the coefficients a have their maximum positive values, i.e., y/S. 
In this case the input function was a constant of amplitude y/S and since the filter has unit gain for D.C., the 
output is the same. Hence the output ensemble has a peak power S. 
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The entropy of the output ensemble can be calculated from that of the input ensemble by using the 
theorem dealing with such a situation. The output entropy is equal to the input entropy plus the geometrical 
mean gain of the filter: 

f w <> r w /W-f\2 

lo ^ C2d f = l 0 ^{-W L ) d f=- 2W - 

Hence the output entropy is 

Wlog4S-2W = Wlog^ 

e L 

and the channel capacity is greater than 

2 S 
Wlog-j 

s 

We now wish to show that, for small — (peak signal power over average white noise power), the channel 
capacity is approximately 

C=Wlog(l + !). " 

More precisely cj Wlog^l 4- -> 1 as ^ -+ 0. Since the average signal power P is less than or equal 

S 

to the peak 5, it follows that for all — 

N 



C<Wlog(l + 0 <Wlog^l + j0. 



Therefore, if we can find an ensemble of functions such that they correspond to a rate nearly W log 



K) 



and are limited to band W and peak S the result will be proved. Consider the ensemble of functions of the 
following type. A series of t samples have the same value, either or -y/S, then the next t samples have 
the same value, etc. The value for a series is chosen at random, probability \ for +\/S and \ for -y/S. If 
this ensemble be passed through a filter with triangular gain characteristic (unit gain at D.C.), the output is 
peak limited to ±S. Furthermore the average power is nearly 5 and can be made to approach this by taking t 
sufficiently large. The entropy of the sum of this and the thermal noise can be found by applying the theorem 
on the sum of a noise and a small signal. This theorem will apply if 

is sufficiently small. This can be ensured by taking — small enough (after t is chosen). The entropy power 
will be S + N to as close an approximation as desired, and hence the rate of transmission as near as we wish 
to 

Wlog 



PART V: THE RATE FOR A CONTINUOUS SOURCE 

27. Fidelity Evaluation Functions 

In the case of a discrete source of information we were able to determine a definite rate of generating 
information, namely the entropy of the underlying stochastic process. With a continuous source the situation 
is considerably more involved. In the first place a continuously variable quantity can assume an infinite 
number of values and requires, therefore, an infinite number of binary digits for exact specification. This 
means that to transmit the output of a continuous source with exact recovery at the receiving point requires, 
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in general, a channel of infinite capacity (in bits per second). Since, ordinarily, channels have a certain 
amount of noise, and therefore a finite capacity, exact transmission is impossible. 

This, however, evades the real issue. Practically, we are not interested in exact transmission when we 
have a continuous source, but only in transmission to within a certain tolerance. The question is, can we 
assign a definite rate to a continuous source when we require only a certain fidelity of recovery, measured in 
a suitable way. Of course, as the fidelity requirements are increased the rate will increase. It will be shown 
that we can, in very general cases, define such a rate, having the property that it is possible, by properly 
encoding the information, to transmit it over a channel whose capacity is equal to the rate in question, and 
satisfy the fidelity requirements. A channel of smaller capacity is insufficient. 

It is first necessary to give a general mathematical formulation of the idea of fidelity of transmission. 
Consider the set of messages of a long duration, say T seconds. The source is described by giving the 
probability density, in the associated space, that the source will select the message in question P(x). A given 
communication system is described (from the external point of view) by giving the conditional probability 
P x (y) that if message x is produced by the source the recovered message at the receiving point will be y. The 
system as a whole (including source and transmission system) is described by the probability function P(x,y) 
of having message x and final output y. If this function is known, the complete characteristics of the system 
from the point of view of fidelity are known. Any evaluation of fidelity must correspond mathematically 
to an operation applied to P(x,y). This operation must at least have the properties of a simple ordering of 
systems; i.e., it must be possible to say of two systems represented by Pi (x,y) and Pz(x,y) that, according to 
our fidelity criterion, either (1) the first has higher fidelity, (2) the second has higher fidelity, or (3) they have 
equal fidelity. This means that a criterion of fidelity can be represented by a numerically valued function: 



whose argument ranges over possible probability functions P{x,y). 

We will now show that under very general and reasonable assumptions the function v(P(x,y)) can be 
written in a seemingly much more specialized form, namely as an average of a function p{x,y) over the set 
of possible values of x and y: 



To obtain this we need only assume (1) that the source and system are ergodic so that a very long sample 
will be, with probability nearly 1, typical of the ensemble, and (2) that the evaluation is "reasonable" in the 
sense that it is possible, by observing a typical input and output x\ and yi , to form a tentative evaluation 
on the basis of these samples; and if these samples are increased in duration the tentative evaluation will, 
with probability 1, approach the exact evaluation based on a full knowledge of >°(x,y). Let the tentative 
evaluation be p(x,y). Then the function p(x,y) approaches (as T -> «>) a constant for almost all (x,y) which 
are in the high probability region corresponding to the system: 



This establishes the desired result. 

The function p(x,y) has the general nature of a "distance" between x and y. 9 It measures how undesirable 
it is (according to our fidelity criterion) to receive y when x is transmitted. The general result given above 
can be restated as follows: Any reasonable evaluation can be represented as an average of a distance function 
over the set of messages and recovered messages x and y weighted according to the probability P(x,y) of 
getting the pair in question, provided the duration T of the messages be taken sufficiently large. 

The following are simple examples of evaluation functions: 

9 It is not a "metric" in the strict sense, however, since in general it does not satisfy either p(x, y) - p{y } x) or p(x,y) +p(y,z) > p(x,z). 



v(P(x,y)) 




p(x,y) -> v(P(x,y)) 



and we may also write 




since 
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1. R.M.S. criterion. 



v= WO -**))• 

In this very commonly used measure of fidelity the distance function p(x,y) is (apart from a constant 
factor) the square of the ordinary Euclidean distance between the points x and y in the associated 
function space. 

p{x,y) = \[{x{t)-y{t)fdt. 

Frequency weighted R.M.S. criterion. More generally one can apply different weights to the different 
frequency components before using an R.M.S. measure of fidelity. This is equivalent to passing the 
difference x(t) — y(t) through a shaping filter and then determining the average power in the output. 
Thus let 

e(t)=x(t)-y(t) 

and „ 

f{t)= T e{r)k{t-r)dr 

J — oo 

then 

p{x,y) = \[f{t?dt. 



3. Absolute error criterion. 



P(x,y) = j£\x(t)-y(t)\dt. 



4. The structure of the ear and brain determine implicitly an evaluation, or rather a number of evaluations, 
• appropriate in the case of speech or music transmission. There is, for example, an "intelligibility" 

criterion in which p{x,y) is equal to the relative frequency of incorrectly interpreted words when 
message jc(r) is received as y(t). Although we cannot give an explicit representation of p(x,y) in these 
cases it could, in principle, be determined by sufficient experimentation. Some of its properties follow 
from well-known experimental results in hearing, e.g., the ear is relatively insensitive to phase and the 
sensitivity to amplitude and frequency is roughly logarithmic. 

5. The discrete case can be considered as a specialization in which we have tacitly assumed an evaluation 
. based on the frequency of errors. The function p(*,y) is then defined as the number of symbols in the 

sequence y differing from the corresponding symbols in jc divided by the total number of symbols in 

x. 

28. The Rate for a Source Relative to a Fidelity Evaluation 

We are now in a position to define a rate of generating information for a continuous source. We are given 
P(x) for the source and an evaluation v determined by a distance function p(x,y) which will be assumed 
continuous in both x and y. With a particular system P(x,y) the quality is measured by 

V = JJ p(x,y)P(x,y)dxdy. 

Furthermore the rate of flow of binary digits corresponding to P(x,y) is 

We define the rate R\ of generating information for a given quality vi of reproduction to be the minimum of 
R when we keep v fixed at vi and vary P x (y). That is: 

P(x,y) 



Rl =MnJJp( X ,y)lo gl $Ld X dy 
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subject to the constraint: 

vi = J J P{x,y)p{x,y)dxdy. 

This means that we consider, in effect, all the communication systems that might be used and that 
transmit with the required fidelity. The rate of transmission in bits per second is calculated for each one 
and we choose that having the least rate. This latter rate is the rate we assign the source for the fidelity in 
question. 

The justification of this definition lies in the following result: 

Theorem 21: If a source has a rate R\ for a valuation v\ it is possible to encode the output of the source 
and transmit it over a channel of capacity C with fidelity as near v\ as desired provided R\ < C. This is not 
possible if R\ > C. 

The last statement in the theorem follows immediately from the definition of R\ and previous results. If 
it were not true we could transmit more than C bits per second over a channel of capacity C. The first part 
of the theorem is proved by a method analogous to that used for Theorem 1 1. We may, in the first place, 
divide the (x,y) space into a large number of small cells and represent the situation as a discrete case. This 
will not change the evaluation function by more than an arbitrarily small amount (when the cells are very 
small) because of the continuity assumed for p(x,y). Suppose that P\(x,y) is the particular system which 
minimizes the rate and gives R\ . We choose from the high probability /s a set at random containing 

2 (R l +c)T 

members where e -> 0 as T -» ©o. With large T each chosen point will be connected by a high probability 
line (as in Fig. 10) to a set of x's. A calculation similar to that used in proving Theorem 1 1 shows that with 
large T almost all x y s are covered by the fans from the chosen y points for almost all choices of the y 's. The 
communication system to be used operates as follows: The selected points are assigned binary numbers. 
When a message x is originated it will (with probability approaching 1 as T — > <*>) lie within at least one 
of the fans. The corresponding binary number is transmitted (or one of them chosen arbitrarily if there are 
several) over the channel by suitable coding means to give a small probability of error. Since R\ <C this is 
possible. At the receiving point the corresponding y is reconstructed and used as the recovered message. 

The evaluation v\ for this system can be made arbitrarily close to v\ by taking T sufficiently large. 
This is due to the fact that for each long sample of message x(t) and recovered message y(t) the evaluation 
approaches v\ (with probability 1). 

It is interesting to note that, in this system, the noise in the recovered message is actually produced by a 
kind of general quantizing at the transmitter and not produced by the noise in the channel. It is more or less 
analogous to the quantizing noise in PCM. 

29. The Calculation of Rates 
The definition of the rate is similar in many respects to the definition of channel capacity. In the former 

R = Mmf[ P(x,y)\og f}*'*). dxdy 
p x (y)JJ y) B P(x)P(y) 

with P(x) and vi = jj P(x,y)p(x,y)dxdy fixed. In the latter 

with P x (y) fixed and possibly one or more other constraints (e.g., an average power limitation) of the form 
K = fJP(x,y)\(x iy )dxdy. 

A partial solution of the general maximizing problem for determining the rate of a source can be given. 
Using Lagrange's method we consider 



n *' y) Iog pWU) + ^M**) + "WW) 



dxdy. 



50 



The variational equation (when we take the first variation on P(x,y)) leads to 

where A is determined to give the required fidelity and B(x) is chosen to satisfy 

jB(x)e' Xf, ^ ) dx=l. . 

This shows that, with best encoding, the conditional probability of a certain cause for various received 
v, P y (x) will decline exponentially with the distance function p(*,y) between the x and y in question. 

In the special case where the distance function p(x y y) depends only on the (vector) difference between x 
and y, 

p{x,y) = p(*-y) 

we have 

[B(x)e- X/>{x - y) dx=l. 

Hence B(x) is constant, say a, and 

P y (x) = ae~ Xfi(x - y) : 

Unfortunately these formal solutions are difficult to evaluate in particular cases and seem to be of little value. 
In fact, the actual calculation of rates has been carried out in only a few very simple cases. 

If the distance function p(x,y) is the mean square discrepancy between x and y and the message ensemble ; 
is white noise, the rate can be determined. In that case we have 

R = Min [H(x) - H y (xj] = H(x) - MaxH y (x) 

with N = (jc-y) 2 . But the Max H y (x) occurs when y - x is a white noise, and is equal to W\ \ogl-neN where 
W\ is the bandwidth of the message ensemble. Therefore 

'R = W\log2ireQ-W\log2TteN 

= WYiog'j? 

where Q is the average message power. This proves the following: 

Theorem 22: The rate for a white noise source of power Q and band W\ relative to an R.M.S. measure 
of fidelity is 

K = Wilog^ 

where N is the allowed mean square error between original and recovered messages. 

More generally with any message source we can obtain inequalities bounding the rate relative to a mean 
square error criterion. 

Theorem 23: The rate for any source of band W\ is bounded by 

Wilog^<*<Wilog§ 
N N 

where Q is the average power of the source, Q\ its entropy power and N the allowed mean square error. 

The lower bound follows from the fact that the MaxH y (x) for a given {x-y) 2 = N occurs in the white 
noise case. The upper bound results if we place points (used in the proof of Theorem 21) not in the best way 
but at random in a sphere of radius y/Q-N. 
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APPENDIX 5 

Let S\ be any measurable subset of the g ensemble, and S 2 the subset of the / ensemble which gives S\ 
under the operation T. Then 

5i = ts 2 . 

Let H x be the operator which shifts all functions in a set by the time A. Then 

H x Si=H x TS 2 = TH x S 2 

since T is invariant and therefore commutes with H x . Hence if m[S] is the probability measure of the set S 

m[// A S,] = m[TH x S 2 ] = m[H x S 2 ] 
= m[S 2 ] = m[Si] 

where the second equality is by definition of measure in the g space, the third since the / ensemble is 
stationary, and the last by definition of g measure again. 

To prove that the ergodic property is preserved under invariant operations, let S\ be a subset of the g 
ensemble which is invariant under //\ and let S 2 be the set of all functions / which transform into S\ . Then 

H X S X = H X TS 2 = TH X S 2 = Si 

so that H X S 2 is included in 52 for all A. Now, since 

m[H x S 2 ] = m[S\] 

this implies 

H X S 2 = S 2 

for all A with m[S 2 ] / 0, 1. This contradiction shows that S\ does not exist. 

APPENDIX 6 

The upper bound, Ni<N\+N 2 , is due to the fact that the maximum possible entropy for a power N\+N 2 
occurs when we have a white noise of this power. In this case the entropy power is N\+N 2 . 

To obtain the lower bound, suppose we have two distributions in n dimensions p(x,-) and_g(*/) with 
entropy powers N\ and N 2 . What form should p and q have to minimize the entropy power W 3 of their 
convolution r(x,): 

r(xi)= J p(yi)q{xi-yi)dy h 

The entropy H3 of r is given by 

#3 = -J r{xi)\ogr(xi)dxi. 
We wish to minimize this subject to the constraints 

Hi = -J p(xi)\ogp{xi)dxi 
#2 = - / q(xi)\ogq{xi)dxi. 



52 



We consider then 

U = -f [r(x)logr(x) + Ap(x)logp(x) + fiq{x)logq(x)] dx 
5U = -J[[1 + \ogr(x)]5r(x) + A[l + logp(x)]*p (x) + + log«(x)] j 9 (x)] dx. 

If /?(x) is varied at a particular argument x,- = the variation in r(x) is 

6r(x) = g(x,- - Si) 

and 

51/ = - J q{xi - Si) logr(x,-) dx; - A log /?(*/) = 0 
and similarly when q is varied. Hence the conditions for a minimum are 

J q(xi-Si)\ogr(xi)dxi = -Alogp(^) 

J p( Xi -Si)\ogr(xi)dxi = -fj,\ogq(si). 

If we multiply the first by p(si) and the second by q(st) and integrate with respect to si we obtain 

H 3 = -\H l 

or solving for A and \i and replacing in the equations 

#1 j q{xi - Si) log r(xi) dx t = -H 3 logp(s f ) 

#2 J P(xi - Si) log r(xi) dxi = -tf 3 log <?(.*,■) . 
Now suppose p(xi) and q(xi) are normal 

U--K 2 
\B-\ n / 2 

q(Xi) = ^^ eXP ^ BijXiXj ' 
Then r(x,) will also be normal with quadratic form Qj. If the inverses of these forms are ay, bij, cy then 

Cij = aij + bij. 

We wish to show that these functions satisfy the minimizing conditions if and only if aij = Kb l} and thus 
give the minimum H3 under the constraints. First we have 

n 1 

logr(x,) = - log — \Qj\ - \ YjQjXiXj 

/n 1 
q{xi - Si) \ogr(xi)dxi = - log — \Qj\ - j ^QjSiSj - \ ^Qjbij. 

This should equal 



#3 
Hi 



which requires Ay = jj-Qj. In this case Aij = 77- fly and both equations reduce to identities. 



« 3 " ' « 2 
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APPENDIX 7 



The following will indicate a more general and more rigorous approach to the central definitions of commu- 
nication theory. Consider a probability measure space whose elements are ordered pairs (x,y). The variables 
x, y are to be identified as the possible transmitted and received signals of some long duration T. Let us call 
the set of all points whose x belongs to a subset S\ of x points the strip over S\ , and similarly the set whose 
y belong to 52 the strip over 52. We divide jc and y into a collection of non-overlapping measurable subsets 
Xj and Yi approximate to the rate of transmission R by 

where 

P(Xi) is the probability measure of the strip over X, 

P(Yi) is the probability measure of the strip over Yi 

P(Xi, Yi) is the probability measure of the intersection of the strips. 

A further subdivision can never decrease R\ . For let X\ be divided into X\ = X[ + X" and let 

P(Y\)=a P{X\)=b + c 

P(X[)=b P(X[,Yx)=d 
P(X[')=c P(Xf,Ki) = e 

P{X u Yi)=d + e. 

Then in the sum we have replaced (for the X\ , Y\ intersection) 

d + e d e 

{d + e)\og 77 ■ : by d\og— -helog — . 
v 7 a{b + c) ab ac 

It is easily shown that with the limitation we have on b, c y d, e, 

'd + e 



b + c 



d+e d d e e 



and consequently the sum is increased. Thus the various possible subdivisions form a directed set, with 
R monotonic increasing with refinement of the subdivision. We may define R unambiguously as the least 
upper bound for R\ and write it 



P(x)P(y) 

This integral, understood in the above sense, includes both the continuous and discrete cases and of course 
many others which cannot be represented in either form. It is trivial in this formulation that if x and u are 
in one-to-one correspondence, the rate from u to y is equal to that from x to y. If v is any function of y (not 
necessarily with an inverse) then the rate from x to y is greater than or equal to that from x to v since, in 
the calculation of the approximations, the subdivisions of y are essentially a finer subdivision of those for 
v. More generally if y and v are related not functionally but statistically, i.e., we have a probability measure 
space (y,v), then R{x, v) < R(x,y). This means that any operation applied to the received signal, even though 
it involves statistical elements, does not increase R. 

Another notion which should be defined precisely in an abstract formulation of the theory is that of 
"dimension rate," that is the average number of dimensions required per second to specify a member of 
an ensemble. In the band limited case 2W numbers per second are sufficient. A general definition can be 
framed as follows. Let f a (t) be an ensemble of functions and let prLMOi //*(')] be a metric measuring 
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the "distance" from f a to fp over the time T (for example the R.M.S. discrepancy over this interval.) Let 
N(e,$,7) be the least number of elements / which can be chosen such that all elements of the ensemble 
apart from a set of measure 5 are within the distance e of at least one of those chosen. Thus we are covering 
the space to within c apart from a set of small measure 6. We define the dimension rate A for the ensemble 
by the triple limit 

A = Li m LimLi,n ^Mll. 

This is a generalization of the measure type definitions of dimension in topology, and agrees with the intu- 
itive dimension rate for simple ensembles where the desired result is obvious. 
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